Mon. Not. R. Astron. Soc. 000,[T]-?? (2009) Printed 12 February 2013 (MN OTeX style file v2.2) 



ATLASGAL — Environments of 6.7 GHz methanol masers 

J. S. Urquhart^*, T. J. T. Moored F. Schuller^ F. Wyrowski\ K. M. Menten\ 
M. A. Thompson'^, T. Csengeri\ C. M. Walmsley^'^, L. Bronfman^, C. Konig^ 

' Max-Planck-Institut fur Radioastronomie, Aufdem Hiigel 69, Bonn, Germany 

^Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Egerton Wharf, Birkenhead, CH41 ILD, UK 
^European Southern Obserx'atory, Alonso de Cordova 3107, Vitacura, Santiago, Chile 

"^Science and Technology Research Institute, University of Hertfordshire, College Lane, Hatfield, ALIO 9AB, UK 

^ Osservatorio Astrofisico di Arcetri, Largo E. Fermi, 5, 50125 Firenze, Italy 

^Dublin Institute for Advanced Studies, Burlington Road 10, Dublin 4, Ireland 

^ Departamento de Astronomi'a, Universidad de Chile, Casilla 36-D, Santiago, Chile 

Accepted ??. Received ??; in oiiginal fonn ?? 



ABSTRACT 

Using the 870-yum APEX Telescope Large Area Survey of the Galaxy (ATLASGAL), we 
have identified 577 submillimetre continuum sources with masers from the methanol multi- 
beam (MMB) survey in the region 280° < £ < 20°; \ b\ < 1.5°. 94 percent of methanol 
masers in the region are associated with sub-millimetre dust emission. We estimate masses 
for ~450 maser-associated sources and find that methanol masers are preferentially associ- 
ated with massive clumps. These clumps are centrally condensed, with envelope structures 
that appear to be scale-free, the mean maser position being offset from the peak column den- 
sity by + 4". Assuming a Kroupa initial mass function and a star-formation efficiency of 
~30 per cent, we find that over two thirds of the clumps are likely to form clusters with masses 
>20Mq. Furthermore, almost all clumps satisfy the empirical mass-size criterion for massive 
star formation. Bolometric luminosities taken from the literature for ~100 clumps range be- 
tween ~100 and 10'' L©. This confirms the link between methanol masers and massive young 
stars for 90 per cent of our sample. The Galactic distribution of sources suggests that the star- 
formation efficiency is significantly reduced in the Galactic-centre region, compared to the 
rest of the survey area, where it is broadly constant, and shows a significant drop in the mas- 
sive star-formation rate density in the outer Galaxy. We find no enhancement in source counts 
towards the southern Scutum-Centaurus arm tangent at f ~ 315°, which suggests that this arm 
is not actively forming stars. 

Key words: Stars: formation - Stars: early-type - Galaxy: structure - ISM: molecules - ISM: 
submillimetre. 



1 INTRODUCTION 

Massive stars (> 8Mq and >10^ Lq) play a hugely important role 
in many astrophysical processes fro m the formation of the first 
solid material in the early Universe dOunne et al.ll2003l) : to their 
substantial influence upon the evoluti on of their host g alaxies and 
future generations of star formation dKennicutdlioOsI) . Given the 
profound impact massive stars have, not only on their local envi- 
ronment, but also on a Galactic scale, it is crucial to understand the 
environmental conditions and processes involved in their birth and 
the earliest stages of their evolution. However, massive stars form 
in clusters and are generally located at greater distances than re- 
gions of low-mass star formation and therefore understanding how 
these objects form is observationally much more challenging (see 
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IZinnecker & Yorkel2007l for a review). Moreover, massive stars are 
rare and they evolve much more quickly than low-mass stars, reach- 
ing the main sequence while still deeply embedded in their natal 
environment. As a consequence of their rarity and relatively short 
evolution large spatial volumes need to be searched in order to iden- 
tify a sufficient number of sources in each evolutionary stage. Only 
then can we begin to understand the processes involved in the for- 
mation and earliest stages of massive star formation. 

There have been a number of studies over the last two 
decades or so that have been used to identify samples of em- 
bedded youn g massive st ars utilizing the presence of methanol 
(e.g.. IWalsh et al.l [19971). IRAS . M SX or GLIMPSE in- 



frared colours (e.g., Molinari et al.l 19961 an d Bronfman et aljl9^3 , 



compact radio emission (e.g., 



iLumsden et aljfaoa and iRobitaille et alj 20081. respectively ) and 

IWood & Churchwelll 119891) 



Al- 



though these methods have had some success, they tend to focus 
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on a particular evolutionary type and in the case of the infrared 
colour selected samples are biased away from complex regions as 
a result of confusion in the images due to the limited angular res- 
olution of the surveys. This is particularly acute for IRAS selected 
samples. Some of these surveys suffer from significant biases. For 
example, methanol masers require a strong mid-infrared source for 
the creation of high enough methanol abundances and to pump the 
maser transitions, while a hot ionising star must be present to ionize 
an UC H n region. Surveys such as these preclude the possibility of 
identifying the very earliest pre-stellar clumps that would need to 
be included in any complete evolutionary sequence for massive star 
formation. 

All of the earliest stages of massive star formation take place 
within massive clumps of dust and gas, which can be traced by 
their thermal dust continuum emission. Dust emission is generally 
optically thin at (sub)millimetre wavelengths and is therefore an 
excellent tracer of column density and total clump mass. As well as 
including all of the embedded stages, dust emission observations 
are also sensitive to the colder pre-stellar phases and so provide a 
means to study the whole evolutionary sequence of the massive star 
formation. Until recently there were no systematic and unbiased 
surveys of dust emission. Most studies to date that have been under- 
taken consiste d of targeted observations of IRAS or maser selected 
samp l es (e.g..lSridharan et al.ll2002l : iFaundez et al.l2004 iHill et alj 



I2OO5I : [Thompson et alj|2006h . and so suff^er from the same prob- 
lems mentioned in the previous paragraph, or have concentrated on 
a single, often exceptionally rich region (e.g., .Motte et al.i,2007l '). 
which may not be representative. However, there are now two large 
unbiased (sub)millimetre surveys available: the AP EX Telescope 
Large Area Survey of the Galaxy (ATLASGAL; ISchuller etall 



20091) at 870 um an d the Bolocam Galactic Plane Survey (BGPS; 
Aguirre et al.l201 iT) at I . I mm. These surveys have identified many 



thousands of sources across the Galaxy with which to compile the 
large samples of massive clumps required to build up a comprehen- 
sive understanding of massive star formation, and test the predic- 



itive accretion 1 


Bonnell et alj|l99ll200ll) 


dMcKee & TarJ 


20031)). 



The complete ATLASGAL ca talogue consists of a pproxi- 
mately 12,000 compact sources (see IContreras et aLll2013l for de- 
tails) distributed across the inner Galaxy. This is the first of a series 
of papers that will investigate the dust properties and Galactic dis- 
tribution of massive star formation. Here we use the association of 
methanol masers, that are considered to be an excellent trace r of 
the early stages of massive star formation iMinier et al.ll2003l) . to 
identify a large sample of massive clumps. 

The structure of this paper is as follows: in Sect. 2 we provide 
a brief summary of the two surveys used to select the sample of 
massive star forming clumps. In Sect. 3 we describe the matching 
procedure and discuss sample statistics, while in Sect. 4 we derive 
the physical properties of the clumps and their associated masers. 
In Sect. 5 we evaluate the potential of the clumps to fom massive 
stars and investigate their Galactic distribution with reference to 
the large scale structural features of the Milky Way. We present a 
summary of the results and highlight our main findings in Sect. 6. 



2 SURVEY DESCRIPTIONS 
2.1 ATLASGAL Survey 

The A PEX Telescope Larg e Area Survey of the Galaxy (ATLAS- 
GAL; ISchuller et al.l |2009|) is the first systematic survey of the 



inner Galactic plane in the submillimeter wavelength range. The 
survey was carried out with the Large APEX Bolometer Camera 
(LABOCA: [siringo et al.ll2009l) . an array of 295 bolometers ob- 
serving at 870 /jm (345 GHz). The 12-metre diameter telescope 
affords an angular resolution of I9"2 FWHM. The initial survey 
region covered a Galactic longitude region of 300° < ( < 60° and 
\b\ < 1.5°, but this was extended to include 280° < I < 300°, 
however, the latitude range was shifted to -2° < b < i° to take 
account of the Galactic warp in this region of the plane and is not 
as sensitive as for the iimer Galaxy (~60 and 100 mjy beam"' for 
the 300° < ^ < 60° and 280 ° <( < 300° regions, respectively). 

IContreras et al.l (1201 3h produced a compact source catalogue 
for the central part of the survey region (i.e., 330 ° < { < 21°) using 
the s ource extraction algorithm SExtractor jBertin & AmoutsI 
fl99^ . Signal-to-noise maps that had been filtered to remove the 
large scale variations due to extended diffuse emission were used 
by SExtractor to detect sources above a threshold of 3cr, where 
cr corresponds to the background noise in the maps. Source param- 
eters were determined for each detection from the dust emission 
maps. This catalogue consists of 6,639 sources and is 99 per cent 
complete at ~6it, which corresponds to a flux sensitivity of 0.3- 
0.4 Jy beam"'. We have used the same source ex traction algorithm 
and method described bv lContreras et al.l ( l20I3h to produce a cata- 
logue for the currently unpublished 280° < { < 330° and 2\° < £ < 
60° regions of the survey. When the source s identified in these re- 
gions are combined with those identified bv lContreras et al.l ( l20I3l) 
we obtain a final compact source catalogue of some 12,000 sources 
(full catalogue will be presented in Csengeri et al. 2013 in prep.). 
The telescope has an r.m.s. pointing accuracy of ~2", which we 
adopted as the positional accuracy for the catalogue. This catalogue 
provides a complete census of dense dust clumps located in the in- 
ner Galaxy and includes all potential massive star forming clumps 
with masses greater than 1,000 Mg out to 20kpc. 

2.2 MMB Survey 

Methanol masers are well-known indicators of the early phases of 
high-mass star formation, in particular source s showing em ission 
in the strong 6.7 GHz Class-II maser transition ( lMentenll99 ih . The 
Methanol Multibeam (MMB) Survey mapped the Galactic plane 
for this maser transition using a 7-beam receiver on the Parkes 
telescope with a s ensitivity of 0.17Jybeam~' and a half-power 
beamwidth of 3.2' l lGreenll2009h . All of these initial maser detec- 
tions were followed up at high-resolution (~2") with the Australia 
Telescope Compact A rray (ATCA) to obta in sub-arcsec positional 
accuracy (0.4" r.m.s.; ICaswell et alJlIoTol) . To dat e the MMB is 
complete between 186° < £ < 20° and \b\ < 2° 1 Caswell et al] 
2010; Gree n et al.l2O10l;ICaswell et al.l201ll;lGreen et al.l2012b and 
has reported the positions of 707 methanol maser sites; these sites 
consist of groups of maser spots that can be spread up to 1" in size 
but are likely to be associated with a single object. 

The MMB catalogue gives the velocity of the peak compo- 
nent and the flux density as measured from both the high sen- 
sitivity ATCA follow-up observations and those measured from 
the initial lower sensitivity Parkes observations. These observa- 
tions were taken over different epochs up to two years apart and 
consequently the measured values can be affected by variability 
of the maser. We have used the velocity and peak flux densities 
measured from the ATCA data for all MMB sources except for 
MMBG321.704-l-0I.I68 as it was not detected in the ATCA obser- 
vations; for this source we use the values recorded from the Parkes 
observations. 
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Figure 1. Examples of the local mid-infrared environment found towards the ATLASGAL-MMB candidate associations. The left panels show examples of 
rejected matches, while the centre and right panels show genuine associations (see text for details). These images are composed of the GLIMPSE 3.6, 4.5 
and 8.0/im IRAC bands (coloured blue, green and red respectively) which are overlaid with grey contours showing the ATLASGAL 870 ^m emission. The 
positions of the methanol masers in each field are indicated by the green circles and the white hatched circle in the lower left comer of each image shows the 
size of the APEX beam at 870 //m. The contour levels start at 2o" and increase in steps set by a dynamically determined power- law of the form D = 3xN' + 2, 
where D is the dynamic range of the submilhmetre emission map (defined as the peak brightness divided by the local rm.s. noise), N is the number of contours 
used (8 in this case), and i is the contour power-law index. The lowest power-law index used was one, which results linearly spaced contours starting at 2(T 
and increasing in steps of 3(T (seel'Thompson et al. ( 2006) for more details). The advantage of this scheme over a linear scheme is its ability to emphasize both 
emission from diffuse extended structures with low surface brightness and emission from bright compact sources. 



3 ATLASGAL-MMB ASSOCIATIONS 
3.1 Matching statistics 

Of the 707 methanol masers currently identified by the MMB sur- 
vey 671 are located in the Galactic longitude and latitude range 
surveyed at 870 yum by the APEX telescope as part of the ATLAS- 
GAL project. This represents ~95 percent of the entire published 
MMB catalogue. As a first step to identify potential matches be- 
tween the methanol masers and the ATLASGAL sources we used a 
matching radius of 120", which is the maximum radius of so urces 
found in the ATLASGAL catalogue i jContreras e t al.l 1201 3h . We 
used the peak 870 yum flux as the ATLASGAL source position 
for these matches. In cases where a methanol maser was found 
to be located within this search radius of two or more ATLAS- 
GAL sources the nearest submilhmetre source was selected as the 
most likely association. This simple radius search identified 637 
potential ATLASGAL-MMB associations from the possible 671 
methanol masers. To verify that these associations are genuine, we 
extracted 3' x 3' regions from the ATLASGAL emission maps and 
inspected these by eye to confirm the emission region and posi- 



tion of the methanol maser are coincident. We failed to find an AT- 
LASGAL source at the MMB position towards nine MMB sources 
and therefore removed these from the associated sample; these 
masers form part of the sample of unassociated masers discussed 
in Sect.[331 

For the potential associations located in the pa rt of the Galaxy 
surve yed by GLIMPSE (i.e., £ > 295 and \b\ < h lBeniamin et alj 
I2OO3I) we extracted 3' x 3' mid-infrared datasets from the project 
archive and produced false colour images of their environments. 
We present a sample of these images in Fig.[T] overplotted with 
contours of the submilhmetre emission. In the left panels of Fig.[T| 
we present some examples of the rejected matches. The final sam- 
ple consists of 628 methanol masers that are positionally coinci- 
dent with 577 ATLASGAL sources (see middle and right panels 
of Fig.[T] for examples of these associations), with two or more 
methanol masers found toward 44 clumps. This sample includes 
~94 per cent of the MMB masers located in the surveys' overlap re- 
gion, however, this represents only ~7 per cent of the ATLASGAL 
sources in the same region. 

For two of the ATLASGAL sources associated with mul- 
tiple methanol masers we find the maser velocities disagree 
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Figure 2. Two-dimensional distribution of the angular offsets between the peak of the ATLASGAL dust emission and the matched MMB source is shown in the 
left panel. The dashed vertical and horizontal lines indicate the x and v = axes, respectively and the blue cross shows the mean offset in longitude and latitude 
of the whole sample. In the right panel we present a histogram showing the surface density of methanol masers as a function of separation between them and 
the position of peak emission of the associated ATLASGAL source. We have truncated the .v-axis of this plot at 40" as there are only 14 ATLASGAL-MMB 
matches with larger separations and the source density effectively falls to zero. The bin size is 2" . 



by more that would be expected if the masers were associated 
with the same molecular complex (i.e., |Av| > 10km s"'). 
These are: AGAL355.538-00.104 which is associated with 
MMB355.538-00.105 (3.8 km s"') arid MMB355.545-00.103 
(-28.2kms-'); arid AGAL313.766-00.862 which is as- 
sociated with MMB313.767-00.863 (-56.3kms-') and 
MMB3 13.774-00.863 (-44.8 km s"'). For each of these we 
have adopted the velocity of the maser with the smallest angular 
offset from the submillimetre peak. 

In the left panel of Fig.|2] we present a plot showing the off- 
sets in Galactic longitude and latitude between the peak positions 
of the ATLASGAL and methanol maser emission. The positional 
congelation between the two tracers is excellent with the mean (in- 
dicated by the blue cross) centred at zero {M = -0.3" ± 0.4" and 
Ab = +0.1S" ± 0.37"). In the right panel of Fig.|2]we present a plot 
showing the MMB angular surface density as a function of angular 
separation between the maser position and the peak of the 870 yum 
emission of the associated ATLASGAL source. This plot reveals a 
strong correlation between the position of the methanol maser and 
the peak submillimetre continuum emission. The distribution peaks 
at separations less than 2" and falls off rapidly as the angular sep- 
aration increases to ~12" after which the distribution flattens off 
to an almost constant background level close to zero. We find that 
~87 per cent of all ATLASGAL-MMB associations have an angular 
separation < 12" (which corresponds to 3<t, where cr is the standard 
deviation of the offsets weighted by the suiface density). The high 
concentration and small offsets between the masers and the peak 
dust emission reveals that the methanol masers are embedded in the 
brightest emission parts of these submillimetre clumps. This would 
suggest that the protostars giving rise to these methanol masers are 
preferentially found towards the centre of their host clumps. This 
is in broad agr eement with the prediction s of the competitive ac- 
cretion model jBonnell et ai]|l99l l 200 ll) where the deeper gravi- 
tational potential at the centre of the clump is able to significantly 
increase the gas density by funneling material from the whole cloud 
toward the centre. However, high angular resolution interferometric 



observations (e.g., with ALMA) are required to conclusively prove 
this hypothesis by measuring the density distribution for the clumps 
on small spatial scales. 

The larger angular offsets are found for approximately 80 
ATLASGAL-MMB associations, which could indicate the pres- 
ence of clumpy substructure that has not been properly identified 
by SExtractor. Alternatively, these matches with large offsets could 
be the result of a chance alignment of a nearby clump with an MMB 
source that is associated with more distant dust clump that falls be- 
low the ATLASGAL sensitivity limit. In the right panels of Fig.[T] 
we present mid-infrared images towards two sources where the 
methanol maser is offset from the peak position of the dust emission 
by more than 12". Inspection of the dust emission (shown by the 
contours) does reveal the presence of weak localised peaks close 
to that of the methanol maser position, which would suggest that 
these two particular sources do possess substructure that has not 
been identified independently by the ATLASGAL source extrac- 
tion method. These two sources are fairly typical of the matches 
with offsets greater than 12" and therefore we would conclude that 
this is the most likely explanation, however, higher sensitivity ob- 
servations are required to confirm this, and to properly characterise 
the dust properties of these methanol maser sites. 

In Fig.[3]we present a plot comparing the methanol flux den- 
sity and the 870 yum peak flux density of the ATLASGAL-MMB 
associations. There is no apparent correlation between these param- 
eters from a visual inspection of this plot, however, the correlation 
coefficient is 0.19 with a significance value of 5x10"* and so there 
is a weak coixelation. Methanol masers are found to be associated 
with a range of evolutionary stages of massive star formation (i.e., 
from the hot molecular core (HMC) through to the ultra-compact 
(UC) HII region stage). Therefore the low level of correlation could 
simply be a reflection of the spread in evolutionary stages covered 
by the methanol masers. It is also important to bear in mind that 
these methanol masers are likely to be associated with the circum- 
stellar envelope/disk of a single embedded source, and given the 
resolution of the ATLASGAL survey it is almost certain that the 
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Figure 3. Comparison of the methanol maser flux density and the peak 
870 (Um flux density of the associated ATLASGAL source. The colours indi- 
cate whet her or not a particular m aser has been associated with mid-infrared 
enussion iGallawav etaDl2013h : red, purple and black colours correspond 
to infrared bright, infrared dark, and sources for which no data is available, 
respectively. The dashed-dotted blue and green lines indicate the 3<t sensi- 
tivities of the ATLASGAL and MMB surveys, respectively. The results of 
a Pearson correlation test returns a coefficient value of 0. 19. 



measured submillimetre flux is related to the mass of the whole 
clump, which is likely to go on to form an entire cluster. It is there- 
fore not surprising that we find the flux densities of the masers and 
the dust clumps to be only weakly coiTelated. 

To investigate whether there is a correlation between the sub- 
millimetre and methanol maser fluxes and the mid-infrared prop- 
erties of the source we have cr oss-correlated the matched sources 
with the lGallawav et al] l l2013h catalogue. Using this catalogue we 
separate our matched sample into three groups, those associated 
with mid-infrared emission, infrared dark sources and those located 
outside the region covered by the GLIM PSE Legacy project, upon 
which the work o f lGallawavet"ai] ( l2013l) is based. The distribution 
of these three groups are shown in Fig.|3]as red, purple and black 
symbols, respectively. Comparing the flux distributions of the in- 
frared bright and dark samples with a Kolmogorov-Smimov (KS) 
test we do not find them to be significantly different. 



3.2 ATLASGAL 870 fim flux distribution 

In Fig.|4] we present plots of the 870 /jm peak and integrated flux 
distribution (upper and lower panels, respectively) of the ATLAS- 
GAL catalogue (grey histogram) and ATLASGAL-MMB associa- 
tions (blue histogram). It is clear from these plots that the methanol 
masers are preferentially associated with the brighter ATLASGAL 
sources in the sense that the probability of an association with a 
maser approaches 100 per cent for brighter clumps. This is partic- 
ularly evident in the peak flux distribution, which reveals that only 
a relatively small number of submillimetre sources brighter than 
~7 Jybeam"' are not associated with a methanol maser. It is also 
clear from these plots that there is a stronger congelation between 
the brightest peak flux density ATLASGAL sources and the pres- 
ence of a methanol maser than between the integrated flux and the 
presence of a maser. The integrated flux is a property of the whole 
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Figure 4. Flux density distribution for ATLASGAL sources in the over- 
lap region (grey filled histogram) and the ATLASGAL-MMB associated 
sources (blue histogram). In the upper and lower panels we present his- 
tograms of the peak and integrated flux densities measured for each source, 
respectively. 



clump/cloud, whereas the peak flux is more likely to be associ- 
ated with the highest column density and/or warmest regions of the 
clump where star formation is taking place. Furthermore, some of 
the ATLASGAL sources with the highest integrated fluxes are large 
extended sources, that can have relatively low peak flux densities. 

In Fig.|5]we present a plot showing the ratio of ATLASGAL 
sources found to be associated with a methanol maser as a function 
of peak 870 /jm flux density. Although the en'ors in the ratios for 
the higher flux bins are relatively large, due to the smaller num- 
bers of sources they contain, there is still clearly a strong corre- 
lation between bright submillimetre emission and the presence of 
a methanol maser. Given that the association rate of ATLASGAL 
sources with methanol masers increases rapidly with peak flux den- 
sity (~I00 percent for sources above 20Jybeam"') an argument 
can be made for the maser emission being effectively isotropic. Al- 
though the radiation beamed from individual maser spots is highly 
directional, there are many very high resolution studies that have 
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Figure 5. Fraction of ATLASGAL sources associated witli a metlianol 
maser as a function of peak 870 fim flux density. Tiie errors are estimated 
using binomial statistics. 



Figure 6. MMB flux distribution for the whole sample (grey filled his- 
togram) and those not associated with an ATLASGAL source (red his- 
togram). The bin size is 0.2 dex. 



revealed significant n umbers of maser sp ots to be associated with 
a single source fe.g.. lGoddietlD l20T ih. These individual maser 
spots are distributed around the central embedded protostar and are 
therefore spatially distinct. 



3.3 Unassociated MMB Sources 

43 MMB masers could not be matched with an ATLASGAL 
source. These masers have integrated flux densities between 0.43 
and 15.65 Jy, with a mean and median value of 2.2 and 1.5 Jy, re- 
spectively. In Fig.|6]we present a histogram showing the distribu- 
tion of maser fluxes for the whole MMB sample (grey filled his- 
togram) and of the unassociated masers (red histogram). The flux 
densities of the unmatched masers are significantly above the MMB 
survey's sensitivity of 0. 17 Jybeam"', and although they are prin- 
cipally found towards the lower flux end of the distribution, they 
are not the weakest masers detected. For comparison the mean and 
median fluxes for the whole maser sample are 47.8 and 5.1 Jy, re- 
spectively. Using a KS test to compare the distribution of the dust- 
less MMB masers with that of the whole MMB catalogue we are 
able to reject the null hypothesis that these are drawn from the same 
population with greater than three sigma confidence. 

It is widely accepted that methanol masers are almost ex- 
clusively associated with high-mass star forming regions (e.g., 
iMinier et al]|2003l : IPandian et aljboiol) . If this is the case then we 
might expect all of these unmatched sources to be located at the 
far side of the Galaxy where their dust emission falls below the 
ATLASGAL detection sensitivity, however, it is possible that these 
are associated with more evolved stars (e.g., IWalsh et al. 200i). We 
will investigate the nature of these unassociated methanol masers in 
more detail in Sect. 5.4. 



4 PHYSICAL PROPERTIES 

In the previous section we identified two subsamples based on the 
possible combinations of ATLASGAL and MMB associations; 1) 
methanol masers associated with the thermal continuum emission 



from cold dust and likely tracing star formation; 2) apparently dust- 
less methanol masers which could be a combination of masers asso- 
ciated with more distant dust clumps not currently detected. These 
two subsamples consist of approximately 94 per cent and 6 per cent 
of the MMB sources in the overlap region, respectively. 

In this section we will concentrate on the ATLASGAL-MMB 
associations to determine the physical properties of their environ- 
ments; these are given for each clump in Table[T| while in Table|2] 
we summarise the global properties. 

4.1 Distances 

Using the velocity o f the peak maser component and a Ga lactic 
rotation model (e.g., iBrand & BliT3ll993l ; keid etal.ir2009h it is 
possible to estimate a particular source's kinematic distance. How- 
ever, for sources located within the solar circle (i.e., < 8.5 kpc of 
the Galactic centre) there is a two fold degeneracy as the source 
velocity corresponds to two distances equally spaced on either 
side of the tangent position. In a follow -up paper to the MMB 
survey iGreen & McClure-Griffithsl boill) used archival HI data 
taken from the Southern and V LA Galact ic Plane Survey s (SGPS 
jMcClure-Griffiths et alj20"05l) and VGPS jStil et alj200"a) . respec- 
tively) to resolve the distance ambiguities to a large number of 
met hanol masers. 

iGreen & McClure-Griffithsl ( l201ll) examined HI spectra for 
734 methanol masers, of which 204 are located at Galactic lon- 
gitudes between ( = 20° -60° and therefore are not included in 
the curr ent MMB cat alogue and are not considered here. How- 
ever, thel Creen & McC lure-Griffiths (2oTl|) study does include 525 
of the 671 MMB sources located in the overlapping ATLASGAL 
and MMB region. Breaking this down further we find that this in- 
cludes 506 of the 628 MMB sources associated with an ATLAS- 
GAL source as discussed in Sect.[3] and 31 of the 43 methanol 
masers not associated with submillimetre emission mentioned in 
Sect.[33] 

iGreen & McClure-Griffithsl ( l201lh used the Galactic rotation 
model of lReid et alj | |2009|) to determine the kinematic distances to 
their sample of methanol masers. However, this model's assump- 
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Table 1. Derived clump parameters. The columns are as follows: (1) and (2) ATLASGAL and MMB names; (3) angular offset between the peak of the 
submillimetre emission and the masers position; (4) ratio of the semi-major to semi-minor sizes of the ATLASGAL source; (5) ratio of the integrated and peak 
submillimetre emission (F-factor); (6) heliocentric distance; (7) Galactocentric distance; (8) effective physical radius; (9) column density; (10) clump mass 
derived from the integrated 870 /.mi emission assuming a dust temperature of 20 K; (1 1) isotropic methanol maser luminosity. 



ATLASGAL name 


MMB Name" 


Offset 


Aspect 


y-factor 


Distance 


Rgc 


Radius 


Log(N(H2)) 


Log(Masscii,n,p) 


Log(LMMB) 






(") 


Ratio 




(kpc) 


(kpc) 


(pc) 


(cm-2) 


(Mo) 


(Jykpc2) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


AGAL28 1.709-01. 104 


MMB281.710-01.104 


2.3 


1.6 


3.44 


4.2 


8.7 


0.50 


22.97 


2.94 


1.88 


AGAL284.352-00.417 


MMB284.352-00.419 


6.5 


1.5 


10.44 


5.2 


8.8 


1.05 


22.84 


3.49 


2.93 


AGAL284.694-00.359 


MMB284.694-00.361 


6.2 


1.0 


2.00 


6.3 


9.2 




22.41 


2.51 


3.22 


AGAL285.339-00.001 


MMB285. 337-00.002 


8.3 


1.7 


5.29 


5.1 


8.7 


0.53 


22.63 


2.96 


3.56 


AGAL287.372+00.646 


MMB287.371+00.644 


6.8 


1.1 


2.59 


5.2 


8.6 


0.19 


22.67 


2.72 


4.56 


AGAL291. 272-00.714 


MMB291.270-00.719 


19.6 


1.4 


11.25 


1.0 


8.2 


0.36 


23.88 


3.13 


2.00 


AGAL291.272-00.714 


MMB291.274-00.709 


20.5 


1.4 


11.25 


1.0 


8.2 


0.36 


23.88 


3.13 


2.94 


AGAL29 1.579-00.432 


MMB291.579-00.431 


4.3 


1.4 


5.05 


8.1 


9.4 


2.02 


23.59 


4.31 


2.92 


AGAL29 1.579-00.432 


MMB291.582-00.435 


14.8 


1.4 


5.05 


7.7 


9.1 


1.90 


23.59 


4.26 


3.31 


AGAL291.636-00.541 


MMB291. 642-00.546 


28.0 


1.5 


14.72 


7.8 


9.2 


2.94 


23.33 


4.49 


2.37 


AGAL29 1.879-00. 809 


MMB291. 879-00.810 


2.7 


1.2 


2.89 


9.9 


10.4 




22.61 


3.27 


3.22 


AGAL292.074-01.129 


MMB292.074-01.131 


8.1 


2.1 


3.18 


3.2 


7.9 




22.34 


2.05 


2.04 


AGAL293. 828-00.746 


MMB293. 827-00.746 


1.4 


1.1 


2.48 


10.7 


10.7 


0.94 


23.10 


3.75 


3.56 


AGAL293.941-00.874 


MMB293.942-00.874 


3.7 


1.2 


3.10 


11.2 


10.9 


0.86 


22.81 


3.59 


3.83 


AGAL294.336-01.705 


MMB294.337-01.706 


5.6 


1.1 


1.72 


1.0 


8.1 




22.59 


1.03 


0.11 


AGAL294.511-01.622 


MMB294.511-01.621 


3.3 


1.1 


3.60 


1.0 


8.1 


0.14 


23.11 


1.89 


1.98 


AGAL294.976-01.734 


MMB294.977-01.734 


3.1 


1.3 


6.73 


0.2 


8.4 


0.05 


23.20 


1.02 


-1.10 


AGAL294.989-01.719 


MMB294.990-01.719 


2.2 


1.1 


4.87 


1.1 


8.1 


0.17 


23.10 


2.03 


2.19 


AGAL296.893-01.306 


MMB296.893-01.305 


2.4 


1.3 


0.98 


10.0 


9.8 




22.52 


2.71 


3.18 


AGAL297.391-00.634 


MMB297.406-00.622' 


68.2 


1.6 


10.82 


10.7 


10.1 


0.50 


22.38 


3.66 


3.31 


AGAL298. 182-00.786 


MMB298. 177-00.795+ 


37.7 


1.0 


3.60 


10.4 


9.9 


1.49 


23.15 


3.94 


3.56 


AGAL298.224-00.339 


MMB298.2 13-00.343"^ 


41.4 


1.3 


9.89 


11.4 


10.5 


3.89 


23.31 


4.61 


3.33 


AGAL298.263+00.739 


MMB298.262+00.739'^ 


2.5 


1.3 


3.23 


4.0 


7.5 


0.37 


22.92 


2.83 


3.47 


AGAL298.63 1-00.362 


MMB298.632-00.362'^ 


6.0 


1.6 


1.31 


11.9 


10.8 




22.31 


2.77 


3.38 


AGAL298.724-00.086 


MMB298.723-00.086' 


5.0 


1.1 


1.84 


10.6 


9.9 


0.65 


22.94 


3.45 


3.20 


AGAL299.012+00.127 


MMB299.013+00.128'^ 


3.5 


1.9 


5.29 


10.2 


9.6 


1.66 


22.72 


3.66 


4.03 


AGAL300.504-00.176 


MMB300.504-00.176+ 


1.4 


1.3 


4.83 


9.6 


9.0 


2.16 


23.11 


3.95 


3.67 


AGAL300.969+01.146 


MMB300.969+01.148* 


7.0 


1.3 


6.49 


4.3 


7.3 


1.25 


23.41 


3.69 


3.05 


AGAL301. 136-00.226 


MMB 30 1.136-00.226' 


2.3 


1.3 


2.59 


4.3 


7.3 


0.78 


23.90 


3.77 


2.59 


AGAL302.032-00.061 


MMB302.032-00.06r 


1.2 


1.2 


3.94 


4.5 


7.2 


0.70 


23.00 


3.10 


3.43 



" Sources with a superscript have been searched for mid-infrared emission by iGallawav et al.l l l2013l) : t and % indicate infrared bright and infrared dark 
sources, respectively, and * identifies the sources they were unable to classify. 

Notes: Only a small portion of the data is provided here, the full table is available in electronic forni at the CDS via anonymous ftp to cdsarc.u-strasbg.fr 
(130.79.125.5) orvia |htipVAasweb.u-strasbgrfr/cgi-bin/qcat?J/MNRAS/| 



Table 2. Summary of derived parameters. 



Parameter 


Number 


Mean 


Standard Error 


Standard Deviation 


Median 


Min 


Max 


Radius (pc) 


375 


1.27 


0.05 


1.01 


0.97 


0.01 


5.69 


Aspect Ratio 


577 


1.51 


0.02 


0.40 


1.40 


1.01 


3.39 


y-factor 


577 


5.63 


0.16 


3.89 


4.55 


0.98 


24.66 


Log[Clump Mass (Mq)] 


442 


3.27 


0.04 


0.77 


3.36 


-2.00 


5.43 


Log[Column Density (cm^^)] 


577 


22.86 


0.04 


1.06 


22.88 


4.16 


24.74 


Log[LMMB (Jykpc-)] 


442 


3.48 


0.05 


0.99 


3.47 


-1.30 


6.34 



tion of a flat rotation curve with a high rotational velocity leads 
to very noticeable differences in the fourth quadrant of the Galaxy 
between the model-derived ta ngent velocities and the emp i rically 
derived values determined by iMcClure-Griffiths & DickevI | |2007|) 
from the HI termination velocities (see their Fig. 8). This leads to 
the near and far distance being located farther from the tangent po- 
sition than would otherwise be expected, and produces a large hole 
in the Gal actic distribution around the tangent positions (see Fig. 4 
of lGreen & McClure-Griffiths..201 1) . 



In order to avoid this w e have used the Galactic rotation 
curve of iBrand & Blitj I I993h as its model tangent velocities are 
a much closer match to the HI termination velocities. The use of 
a different rotation model in most cases results in only a slight 
change in the estimated kinematic distances from those given by 
iGreen & McClure-Griffithsl jlOIlh . In addition to the minor dif- 
ference in kinematic distances imposed by the change of rotation 
model (~1 kpc) we also: 1) place any source with a velocity within 
10 km s"' of the tangent point at the tangent distance, and 2) place 
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Heliocentric Distance (kpc) 

Figure 7. Heliocentric distance distribution for tlie 473 MMB distances 
available (grey filled histogram) and the distance of 3 1 MMB sources not 
associated with an ATLASGAL source (green histogram). The bin size is 
1 kpc. 



any sources within the solar circle at the near distance if a far- 
distance allocation would lead to a height above the mid-plane 
larger than 4 times the scale heigh t of y oung massive stars (i.e., 
-30 pc; Reed 2000; Urquhart et al. 2011ah . 
[Green & McClure-Griffithsl (|2oTll) 

provides distance solu- 
tions for 385 of the ATLASGAL-MMB associations and 16 of the 
unassociated MMB sources. We have adopted their distances for 
378 of the ATLASGAL-MMB associations and find the values 
determined from the difference rotation curves, after applying 
the two criteria mentioned in the previous paragraph, agree 
within 1 kpc in ever y case. We only disagree with th e distance 
allocations given by [Green & McClure-Griffithsl ( [201 l[) for seven 
sources. Of these, we have associated four sources with the 
G305 complex (AGAL305.36U00.186, AGAL305. 362-^00.151, 
AGAL305.799-00.244 and AGAL305.887-hOO.016), three 
located at the tangent position (AGAL309.384-00.134, 
AGAL311.627-h00.266 and AGAL336.916-00.022). Only 
one source (AGA L35 1 .774-0 0.537) is placed at the near distance 
of 0.4 kpc by the iReid et ak 1 2009i) model but is placed outside 
the solar circle by the [Brand & Blitzl ( Il993h model at a distance of 
17.4 kpc. 

There are 83 ATLAS GAL-MMB associations that have not 
been assigned a distance bv [Green & McClure-Griffithij201lh . We 
have allocated distances to 64 of these. Fourteen have been placed 
at the tangent position, fourteen have been associated with a well- 
known complex (i.e., G305 and W31) and 19 have been placed at 
the near distance. Fifteen of the sources placed at the near distance 
are because a far-distance allocation would place than more then 
120 pc from the Galactic mid-plane. Finally, the velocities of sev- 
enteen sources places them outside the solar circle, and thus, these 
sources do not suffer from the kinematic distance ambiguity prob- 
lem. 

In total we have distances to 442 ATLASGAL-MMB associ- 
ated sources and their distribution is shown in Fig.[7](grey filled 
histogram). The distribution is shown to be bimodal, with peaks at 
3-4 kpc and ll-12kpc. The low number of sources between these 



two peaks is due to an almost total lack of any MMB sources within 
3 kpc of the Galactic centre. 

In addition to the distances we have determined for the 
ATLASGAL-MMB associations, we have obtained distances for 33 
of the 43 MMB masers not associated with thermal dust emission. 
Nineteen of these were drawn from [Green & McClure-Griffithi 
l [2011 ), eleven are found to be located in the outer Galaxy using 
the [Brand&Blitj | [l993l) model, two others are located at the tan- 
gent position and one is placed at the near distance as a far distance 
would place more than 120 pc from the Galactic mid-plane. It is 
clear from the distribution of this sample of sources (green his- 
togram shown in Fig.[7]l that the majority have distances larger than 
9 kpc. In the previous section (i.e.. Sect. 3.3) we suggested that one 
possible explanation for the non-detection of submillimetre dust 
emission from these sources could be that they are located at larger 
distances; this distribution offers some support for that hypothe- 
sis. The median value for the distance for these unassociated MMB 
sources is ~ 13 kpc compared with a median value of ~5kpc for 
the ATLASGAL-MMB associations. A KS test on the two samples 
shows their distance distributions are significantly different and we 
are able to reject the null hypothesis that they are drawn from the 
same population. 



4.2 Sizes and morphology 

Using the kinematic distances discussed in the previous section and 
the effective angular radii derived for the clump s by SExtractor 
we ar e able to estimate their physical sizes (see [Rosolowskv et alj 
[2OIOI for definition of effective radius). The distribution of sizes 
is presented in the upper panel of Fig.[8]and ranges from 0.1 to 
several pc, with a peak at ~I pc. The dashed vertical lines shown 
on this plot at radii of 0.15 and 1.25 pc indicate the boundary be- 
tween core s and clumps, and clu mps and clouds , respectively (as 
adopted bv [Dunham et alj[20Ila[ from Table 1 of [Bergin & Tafallal 
[2007). However, we note that the distribution is continuous with 
no features at these sizes so the definitions are probably somewhat 
arbitrary. 

As also reported bv [Tackenberg et all ( [2012b . from a study of 
ATLASGAL candidate starless clumps we find no correlation be- 
tween angular size and distance, which results in an approximately 
linear correlation between physical sizes and distance. It is impor- 
tant to bear in mind that at distances of a kpc or so we are primar- 
ily resolving structures on the size scale of cores, at intermediate 
distances, clumps, and at greater distances, entire cloud structures. 
This will have an effect on some of the derived parameters such as 
column and volume densities. However, although the sample cov- 
ers a large range of sizes, the majority falls into the size range sug- 
gested for clumps, which are more likely to be in the process of 
forming clusters rather than a single massive star. For simplicity 
we refer to our sample as clumps, with the caveat that it includes a 
large range of physical sizes. 

In the lower panel of Fig.[8]we plot the source radius as a func- 
tion of heliocentric distance. The colours of the symbols used in 
this plot give an indication of the aspect ratio of each source, the 
values of which can be read off from the colour bar to the right 
of the plot. This plot clearly illustrates that at larger distances we 
are probing larger physical structures, however, the aspect ratio of 
this sample of objects does not appear to have a significant dis- 
tance dependence. This may suggest that even at larger distance 
the molecular clouds associated with MMB sources are still single 
structures. 

In the upper panel of Fig.|9| we present the distribution of 
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Figure 8. Upper panel: the distribution of eflfective radius for tlie wliole 
ATLASGAL-MMB associated sample (grey filled histogram). Lower panel: 
the radius and aspect ratio distribution as a function of heliocentric distance. 
The left and right (lower and upper) dashed vertical (horizontal) lines in the 
upper panel (lower panel) indicate the radii separating cores and clumps 
(0.125 pc), and clumps and clouds (1.25 pc), respectively. The solid curved 
line shown in the lower panel shows the physical resolution of the survey 
based on the APEX 19.2" beam at 870 ,um while the colours give an indi- 
cation of the aspect ratio of each sources (see colour bar on the right for 
values). 

the aspect ratios of ATLASGAL-MMB sources (blue histogram), 
which is shown against that of the whole population of compact AT- 
LASGAL sources (grey filled histogram). It is clear from this plot 
that the ATLASGAL-MMB sources have a significantly smaller as- 
pect ratio than the general population, which would suggest they 
are more spherical in structure. The mean (median) values for the 
ATLASGAL-MMB and the whole population are 1.51±0.02 (1.40) 
and 1.69±0.01 (1.55), respectively. A similar median value of 1.3 
was found by [Thompson et alj ^20061) from a programme of tar- 
geted SCUBA submillimetre observations of clumps associated 
with UC H II regions. 

Another way to compare the morphology of the ATLASGAL- 
MMB associations with the general population of ATLASGAL 



Figure 9. In the upper and lower panels we present plots showing the source 
aspect ratio and the y-factor, respectively, of the whole ATLASGAL sam- 
ple (grey histogram) and the ATLASGAL-MMB associated sources (blue 
histogram). The ATLASGAL-MMB distributions have been noiTnalised to 
the peak of the whole ATLASGAL sample. The bin sizes used ai'c 0.1 and 
1 for the aspect ratio and the F-factor, respectively. 



sources is to look at the ratio of their integrated to peak fluxes; 
this is referred to as the Y-factor. This has been used to investi- 
gate the general extent of submillimetre clumps associated with 
high-mass pro t ostella r co res (HMPOs) and UC H ii regions by 
IWilhams etal] ( |2004|) and lThompson et all ( l2006h . respectively. In 
the lower panel of Fig.|9] we plot the K-factor for the whole AT- 
LASGAL compact source population (filled grey histogram) and 
the ATLASGAL-MMB associated sources (blue histogram). The 
y-factors for both distributions peak at similar values, between 2-4, 
however, the ATLASGAL-MMB associated sources have signifi- 
cantly lower overall K-factors (with a median value of 4.8 com- 
pared to 6.7 for the whole ATLASGAL compact source popula- 
tion). A KS test is able to reject the null hypothesis that these are 
drawn from the same population with greater than three sigma con- 
fidence. 

Peaks in the F-factor around 3 are also seen in the 
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Figure 10. Left panel: The isothermal dust mass distribution of the ATLASGAL-MMB associated clumps assuming a dust temperature of 20 K. The vertical 
blue line indicates the completeness limit (see text for details). The bin size is 0.25. Right panel: Dust mass distribution of ATLASGAL-MMB associated 
clumps as a function of heliocentric distance. The colours give an indication of the peak column density of each source (see colour bar on the right for values) 
and the solid black line and the grey filled region indicated the mass sensitivity limit of the survey and its associated uncertainty assuming a dust temperature 
of 20±5 K. 



HMPO and UC H n re gion samples of Iwilliams et all (|2004|) and 
Ixhompson et alj ( l2006l) and appear to have a similar distribution to 
our sample sample. The large range of heli ocentric distances over 
which their HMPO sample is distributed led I Williams et alj ( |2004) 
to suggest that the envelope structures may be scale-free. The AT- 
LASGAL clumps are similarly located over a range of heliocentric 
distances and thus the ATLASGAL beam traces structure on scales 
from a few times 0.1 pc to a several pc. As we observe similar Y- 
factors for a range of spatial resolutions this suggests that the radial 
density distribution follows a similar power law over a range of spa- 
tial scales. Hence the structure of the clump envelopes (at least on 
pc scales) are likely to follow a scale-free power law. The fact that 
all three samples (ATLASGAL-MMB, HMPO and UC H ii regions) 
are likely to cover the whole range of embedded star formation 
and broadly show the s ame K- factor properties wo uld seem to sup- 
port th is. Similarly to William s et al. (2004b and [Thompson et al] 
hood) , we find that a significant amount of the mass associated 
with the ATLASGAL-MMB sources is found in the outer regions 
of the clumps, from which we conclude that this situation does not 
change significantly over the evolution of the embedded stars. 

In summary we have determined that overall the ATLASGAL- 
MMB associated sources are roughly spherical, centrally con- 
densed clumps that appear to have a scale-free envelope with a 
methanol maser coincident with the peak of the submillimetre 
emission. 



4.3 Isothermal clump masses 

In calculating masses and column densities we have assumed that 
all of the measured flux arises from warm dust, however, free- free 
emission from embedded ionised gas and/or molecular line emis- 
sion could make a sign ificant contribution f or broadband bolome- 
ters such as LABOCA. ISchuller et al.l ( |2009|) considered these two 
forms of contaminating emission and concluded that even in the 
most extreme case of the giant HII region associated with W43, and 
the CO (3-2) lines associated with extreme outflows, hot cores and 



photon-dominated regions, the contribution from free-free emis- 
sion and molecular lines is of the order 20 and 15 percent, respec- 
tivel y, and in the major ity of cases will be almost negligible (see 



tivel y, and m the ma|or i 
also jPrabek et al.ll2012l) . 



Assuming the dust is generally optically thin and can be 
characterised by a single temperature we are able to estimate 
the isothermal dust masses for t he AT LASGAL-MMB associated 
sources. Following iHildebrandl ( Il983l) , the total mass in a clump 
is directly proportional to the total flux density integrated over the 



M 



D-SyR 



(1) 



where is the integrated 870 yum flux, D is the heliocentric dis- 
tance to the source, R is the gas-to-dust mass ratio (assumed to be 
100), By is the Planck function for a dust temperature T/j, and /<•,, 
is the dust abso rption coefficient take n as 1.85 cm^ g"' (this value 
was derived by ISchuller et al. I (l2009l) by interpolati ng to 870 ;um 
from Table 1, Col. 5 of lOssenkopf & Hennind jl994h ). 

As reliable dust temperatures are not available for the majority 
of our sample, we make the simplifying assumption that all of the 
clumps have approximately the same temperature and set this to be 
20 K. Single-dish ammonia studies have derived kinetic gas tem- 
peratures for a large of number massive star formation regions that 
cover th e full range of evolu tionary stages. These include methanol 
masers jPandian et ai]|2012l) : 1 . 1 mm thermal dust sources identi- 
fied by the Bolocam Galactic Plane S urvey (BGPS; 



Dunham et alj 



1201 Ibl) : 870 ^m ATLASGAL soui -ces jWienen et al .. 
Red MSX Source Survey (RMS: lUrguhart et al.ll201 



Z012|) . 
Mil 



20121) and the 
For the 



methanol masers mean and median kinetic temperatures of 26 K 



All of these observati onal programmes used either the Green B ank 



Telescope (FWHM -30" lUrquhait et al. 



the E flfelsberg telescope (FWHM -40": 

I20TI and therefore have comparable resolution and 



■I2OI Ibl: iDunh am et al.''201 \\h or 
Pmldian et al.l i2012: Wie nen et al] 



sensitivity. 
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and 23.4 K were reported bv lPandian et alj ( l2012h . and for the mas- 
sive young stellar objects (MYSOs) and UC H ii reg ions have mean 
and me dian kinetic temperatures of 22.1 and 2 1 .4 K lUrquhart et al J 
feOlld) . while Wienen et al. ( 20T2I) report kinetic temperatures of 
~24K for a subs ample of ATLASGAL sources associated with 
methanol masers. iDunham et al. report lower kinetic gas 

temperatures 15.6±5.0K, however, their sample includes a larger 
fraction of starless clumps and so the lower mean temperature is 
expected. 

A kinetic gas temperature of ~25 K would seem to char- 
acterise the clumps that show evidence of star formation, how- 
ever, this temperature is likely to be an upper limit to the clump- 
averaged kinetic temperature as these observations were pointed 
at the peak emission of the clumps and the kinetic temperature is 
likely to be significant l y lower towards the e dges of the clumps 
jZinchenko et al J iggTh . fDunham et al.l j2011bh estimate that using 
the peak kinetic temperature for the whole clump may underes- 
timate the isothermal mass by up to a factor of two. Therefore 
we have chosen to use a value of 20 K, c onsistent with a number 
of similar studies (cf. iMotte et allbool Iffill et alj|2005l) . Given 
that the true clump-averaged kinetic temperatures of the sample 
are likely to range between 15 and 25 K the resulting uncertainties 
in the derived isothermal clump masses of individual sources are 
±43 percent (allowing for an uncertainty in temperature of ±5K 
which is is added in quadrature with the 15 per cent flux measure- 
ment uncertainty). However, this is unlikely to have a significant 
impact on the overall mass distribution or the statistical analysis of 
the masses. We note that 10 per cent of the ATLASGAL-MMB as- 
soci ations are also associated with embedded UCHII regions. How- 
ever, IUrquhartetaLl l l2011bl) found that the presence of an UC H 11 
region only results in an increase in clump-averaged kinetic tem- 
peratures of a few Kelvin. 

In the left panel of Fig.[TO] we present a plot of the isother- 
mal dust mass distribution, while in the right panel we show 
the mass distribution as a function of heliocentric distance. It is 
clear from the right panel of this figure that we are sensitive 
to all ATLASGAL-MMB associated clumps with masses above 
1,000 Mo across the Galaxy, and our statistics should be complete 
above this level. (This completeness limit is indicated on the left 
panel of Fig.[TO]by the vertical blue line.) In this regard it is inter- 
esting to note that the mass distribution peaks at several thousand 
solar masses (see left panel of Fig.[TO]l. which is significantly above 
the completeness limit, and so the drop off in the source counts be- 
tween the completeness limit and the peak is likely to be real. This 
is an important point as it confirms that the methanol masers are 
preferentially asso ciated with very mas sive clumps. 

According to lLada & Ladal ( |2003|) and lMotte et"aLl ( l2003h the 
radius and mass required to form stellar clusters is of the order 0.5- 
I pc and 1 00-1, 000 Mq, respectively. Given the sizes and masses 
of the ATLASGAL-MMB associated clumps it is highly likely that 
the majority are in the process of forming clusters. Assuming a 
star formation efficiency (SEE) of 30 per cent and an in itial mass 
function (IME; Kroupa 2001), Tackenb erg et al. I l l20I2h estimate 
that a clump mass of ~ 1,000 Mq is required to have the potential 
to form at least one 20 M© star, while a clump of ~3,000Mo is 
required to form at least one star more massive than 40 M©. It is 
consistent with the assumption that methanol masers are associated 
mainly with high-mass star formation to find that the majority of 
ATLASGAL-MMB associations (~72 per cent) have masses larger 
than ~ 1,000 Mq, and thus, satisfy the mass requirement for massive 
star formation. 

We also note that approximately a third of the ATLASGAL- 




10" 10' 
Clump Mass (Mj) 

Figure 11. Differential mass distribution. The red line shows a linear fit to 
the data and has a slope of — 2.0±0.1 which is close to the Salpeter (1955) 
value and similar to that found for a number of other studies. The eiTors are 
estimated using Poisson statistics and the bin size used is 0.25 dex. 



MMB associations (~28 per cent) have masses below what is 
thought to be required to form at least one massive star. However, 
all of these sources also tend to be more compact objects (in most 
cases <0.3 pc) and may have a higher star formation efficiency and 
be forming smaller stellar sy stems where the stellar IME does not 
apply (e.g.. lMotteetal.ll2007h . 



4.3.1 Peak column densities 

We estimate column densities from the peak ffux density of the 
clumps using the following equation: 



S,R 



(2) 



where CI is the beam solid angle, /i is the mean molecular weight 
of the interstellar medium, which we assume to be equal to 2.8, 
and niH is the mass of an hydrogen atom, while k,, and R are as 
previously defined. We again assume a dust temperature of 20 K. 

The derived column densities are in the range ~10^^"-'' cm"^, 
peaking at 10^' cm"^. This coiTesponds to a surface density of few 
times 0.1 gcm"^, which is a factor of a few lower than the value 
of 1 g cm"^ predicted to be the lower limit for massiv e star for- 
mation (i.e.. lMcKee & Tanll2003l : lKrumholz et alj|200g|) . However, 
we should not read too much into this as the column densities of 
the more distant sources can be affected by beam dilution. This ef- 
fect is nicely illustrated in the right panel of Eig.[TO]where we use 
colours to show the column densities as a function of distance (see 
colour bar for values). There is clearly a dependence of column 
density on distance as, given the resolution of the survey, we are 
sampling larger scale physical structures as the distance increases, 
which preferentially reduces the column densities of more distant 
sources. It is likely that these more distant sources would frag- 
ment into smaller and dens er core-like structures at higher resolu- 
tion (e.g.. lMotte et alj2007l) . Therefore one should exercise caution 
when drawing conclusions from bea m-averaged quantities su ch as 
the column and volume densities (cf. lTackenberg et al.l2012h . 
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Figure 12. In the left panel we present the methanol maser luminosity distribution of the ATLASGAL-MMB associations. The vertical blue dashed line shows 
the MMB luminosity completeness limit. The bin size used is 0.25 dex. In the right panel we show the source luminosity as a function of heliocentiic distance. 
The solid black line indicated the MMB surveys luminosity sensitivity limit. The colours used to plot each source show the clump-averaged volume density; 
see color bar to the right of the plot for values. 



4.3.2 Clump mass function 

In Fig.[TT] we present the differential mass distribution for the 
ATLASGAL-MMB associated clumps. On this plot the dashed ver- 
tical blue line indicates the completeness limit and the solid red 
line shows the results of a linear least-squares fit to the bins above 
the peak in the mass distribution shown in the left panel of Fig.llOl 
(~3,000Mo). This line provides a reasonable fit to all of the bins 
above the peak mass. It does not fit the two bins just above the 
completeness limit, which may suggest that a second power law 
is required to account for these mass bins. The derived exponent 
{a where AN I AM oc M") of the fit to the high-mass tail for the 
ATLASGAL-MMB associated sources is - 2.0 ±0.1, which is sim - 
ilar to values (-2.0 to -2.3) derived bv |Williams et alj ( |2004 . 



iReid & Wilsoi] ( l2005h. iBeltran et al.| (l2006h from studies of more 



evolved stages and iTackenberg et al. 12012 who targeted a sample 
of starless clumps identified from ATLASGAL data. This similarity 
between the clump mass function over the different stages of mas- 
sive star formation would suggest that the CMF does not change 
significantly as the embedded star formation evolves. 



4.4 Methanol maser luminosity 

In the left panel of Fig.[T2]we present the 6.7 GHz luminosity distri- 
bution of the ATLASGAL-MMB associated masers for which we 
have derived a distance. These "luminosities" have been calculated 
assuming the emission is isotropic using: 



AttD-S, 



(3) 



where D is the heliocentric distance in kpc and 5 ,, is the methanol 
maser peak flux density in Jy. Lmmb thus has units of Jykpc-. 
The distribution shown in the left panal of Fig.[T2] peaks at 
~3,000 Jykpc- , but is skewed to higher luminosities with a mean 
value of 3 X lO'' Jy kpc^. In the right panel of Fig.[T2]we show the 
methanol maser luminosities as a function of heliocentric distance. 
We can see from this plot that the MMB survey is effectively com- 



plete across the Galaxy (~20kpc) to all methanol masers with lu- 
minosities greater than ~ 1,000 Jykpc^. This completeness limit is 
shown by the dashed vertical line plotted on the luminosity distri- 
bution shown in the left panel of this figure. 



4.4.1 Luminosity-volume density correlation 

There has been a number of recent methanol and water-maser and 
dust-clump studies that have reported a trend towards lower mean 
volume densities with increasing maser luminosity, which has been 
interpreted as the result of the evolut i on of the embedded star for- 
mation (i.e., lBreen & Ellingsenll201ll : iBreen et al.l2oTTll2012h . 

To test this trend for our sample of methanol masers, we 
plot the maser luminosity as a function of clump-averaged vol- 
ume density in Fig.[T3]for the whole sample of 374 ATLASGAL- 
MMB associations for which we have sufficient data to derive the 
maser luminosity and volume density^ We more or less replicate 
the correlation coe fficient and least-squares fit gradient reported by 
iBreen et al. U20T l|) using the whole sample. Correlation coefficients 
are -0.43 and -0.46 and gradients -0.78 ± 0.09 and -0.85 ± 0.16 
for the fit presented in Fig.[T3]and Fig. 2 of iBreen et al.l l l201ll) . re- 
spectively. However, from a casual inspection of the mean volume 
densities (see colour bar of the right panel of Fig.[T2) it is clear that 
there is a significant distance dependency. It would appear that the 
poorer sensitivity to lower luminosity masers, and the increase in 
the spatial volume being sampled at larger distances, results in a 
decrease in the mean volume-densities of clumps and an increase 
in maser luminosities. 

To test this further, we have performed a partial Spearman 
coiTelation test of the mean volume-densities and maser lumi- 
nosities to remove their mutual dependence on the distance (e.g.. 



^ Beam-averaged volume densities would be significantly larger However, 
we us e clump-averaged v alues here in order to be consistent with the anal- 



ysis of lBreenetal.l201ll 
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Figure 13. The isotropic methanol-maser luminosity as a function of 
clump-averaged volume density. The red lines show the result of a linear 
least-squares fit to the data. 

lYates et al]|l986l) . of the form tab.c, where 

_ ''AB - ''AC^'BC , 

= [(l-.ic)(l-4c)]"^' ^ ^ 

where A, B and C are the maser luminosity, mean particle density 
and distance respectively, and tab, ''ac and r^c are the Spearman 
rank correlation coefficients for each pair of parameters. The signif- 
icance of the partial rank correlation coefficients is estimated using 
rAB,c[(^ - 3) /(l - yApr)]'^^ assumin g it is distributed as Student's 
t statistic (see lCollins & Mani]l 19981 for more details). 

This returned a partial correlation coefficient value of -0.06 
and a p-value ~0.6 and we are therefore unable to exclude the 
null hypothesis that the sample is drawn from a population where 
rho = 0. We would conclude that there is no intrinsic correlation be- 
tween maser luminosity and clump-averaged volume density. We 
additionally performed tests on distance selected subsamples and 
obtained the same result. The scales traced by the ATLASGAL 
observations are much larger than those of either an HMC or an 
UCHII region, which by definition are < 0.1 pc. We do not see 
evidence for radical changes in the density distribution of the en- 
velope on pc scales with evolutionary state (the distribution of Y 
values being similar). Hence, the density distribution of the enve- 
lope does not evolve appreciably over the relevant timescale. For 
the clump-averaged volume density, the volume density of the en- 
velope dominates over the much smaller core, and so we would not 
expect volume densities derived from single dish observations to 
change appreciably either. 

Our conclusion that there is no correlation between clump 
density and methanol maser lum inosity is supported by a recent 
study bv lCvganowski et"al] | |2013|) who also failed to find a correla- 
tion between water maser luminosities and clum p densities towards 
a sam ple of extended green objects (EGOs; ICyganowski et alj 
I2OO8I) . This is not surprising, since for the masers to work, the den- 
sity must be within a narrow range, which has little to do with the 
density sampled in the ATLASGAL beam. Since the H2 number 
densities form a central part of the evolutionary arguments put for- 
ward by Breen et al. the results presented here cast significant doubt 
on some of their conclusions. 



Figure 14. The relationship between clump mass and luminosity of the em- 
bedded maser. The region in the top right comer of the plot outlined by the 
dashed box shows the part of the parameter space we are complete to both 
methanol maser luminosity and clump mass. The dashed blue lines show 
the results of a linear fit to the data. Using a partial correlation function to 
remove the dependency of both of these parameters on distance^ we obtain 
a coefficient value of 0.37. 

4.5 Maser luminosity-clump mass: completeness 

In Fig.[T4]we present a plot of the maser luminosity as a function of 
clump mass for the 442 ATLASGAL-MMB associations for which 
we have a distance. The upper right region outlined by the grey 
dashed line of this figure indicates the region of parameter space 
that we are complete across the Galaxy to both methanol maser 
luminosity and clump mass; this region contains 280 sources. 

Applying a partial coiTelation function to the clump mass and 
maser luminosit y to remove the depend ence of these two parame- 
ters on distance jCollins & MannI 19981) we obtain a partial correla- 
tion coefficient value of 0.37 with a p-value <k 0.01, suggesting that 
there is a weak correlation between the maser luminosity and clump 
mass. We fit these data with a linear least squared function the result 
of which is shown in Fig.[T4]as a blue dashed line. The parameters 
ofthefitare:Log(LMMB) = (0.857±0.046) M,i„,^p-H(0.692 + 0.155). 
Within the en'ors the fit to the data has a gradient close to one and 
so the relationship between mass and maser luminosity is relatively 
linear. As discussed in Sect. 4.3 the most massive clumps are likely 
to be forming more massive stars (assuming a Kroupa IMF and a 
SFE of 30 per cent) and therefore the weak correlation between the 
clump mass and maser luminosity may be related to this. This sug- 
gests that higher (isotropic) maser luminosity is related to higher 
stellar luminosity in some way, perhaps via the pumping mecha- 
nism or maybe the larger clump/core just provides a longer maser 
amplification column. 



5 DISCUSSION 

There have been a number of studies that have tried to firmly es- 
tablish a connection between the presence of a methanol maser and 
ongoing massive star formation and these have been relatively suc- 
cessful. Most have searched for methanol masers towards low lu- 
minosity protostellar sources and, when no masers were detected. 
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have concluded that methanol masers are exclusiv e ly associated 
with high-mass protostars (e.g.. iMinier et al.ll2003t Isourke et alj 
l2005h . However, most of these surveys have focused on small sam- 
ples, with poorly defined selection criteria, and often use IRAS 
fluxes to determine the luminosity of the embedded source and so 
the luminosities may have been overestimated. So, although these 
studies have set a lower limit to the luminosity of the associated 
protostar, it is unclear whether the results obtained are applicable 
to the whole methanol maser population. 

In this subsection we will draw on the results presented in the 
previous sections to test this hypothesis. The ATLASGAL-MMB 
sample presented in this paper includes 94 per cent of the MMB 
sources in the overlap region of the two surveys, and 90 per cent of 
the entire MMB published catalogue. Moreover, since the number 
of methanol masers in the whole Galaxy is not expected to exceed 
a few thousand, our sample is likely to incoiporate a large fraction 
of the whole Galactic population. Therefore any statistical results 
drawn from this sample will be a fair reflection of the properties of 
the general population. 



5.1 Empirical mass-size relationship for massive star 
formation 

5. 1. 1 Criterion for massive star formation 

In two papers iKaufTmann eTal] ( l2010j |hl) investigated the mass- 
radius relationship of nearby (<500pc) molecular cloud complexes 
(i.e., Ophiuchus, Perseus, Taurus and the Pipe Nebula) and found 
that clouds that were devoid of any high-mass star formation gen- 
erally obeyed the following empirical relationship: 

m('-)< 580Mo(i?eff/pc)'-'' (5) 

where R^g is the effective radius as defined by iRosolowskv et"al] 
( |2010|) FI Comparing this mass-size relation with sam ples of known 
high-mass star fo r ming r egions such as tho se of Beuther et ^ 
' 2002h . iHill et alj (|2005|) . iMotte et alj j2007l) and .Mueller et aU 
2002!) . Kaiiffmann et alj ( l2010bl) ound that all of these regions 
occupied the opposing side of the parameter space (i.e., where 
m(r) > 580 Mo (i^eff/pc)' ". This led them to suggest that the 



relation may approximate a requirement for massive star forma- 
tion, where only more massive clumps have the potential to form 
massive stars. However. iKauffmann et alj ( l2010bl) states that larger 
samples of massive star forming clouds are required to strengthen 
this hypothesis. 

In Fig.[T5] we present the mass-size relationship for the 
ATLASGAL-MMB associated sources. This sample consists of 
375 clumps that have distance estimates and are spatially re- 
solved in the APEX beam. Of these, we find that 363 have masses 
larger than the l imiting mass for their size, as determined by 
lKauff"mann et alj ( l2010bl) for massive star formation. This corre- 
sponds to ~97 per cent of the sample and, although this in itself 
does not confirm that the embedded source is a high-mass proto- 
star, it does at least suggest that these clumps have the potential to 
form one. This also supports our earlier statement (in Sect. 4.3) that 



^ Note that when derivin g this relationship | K auffman n et alj J2010bl) re- 
duced the dust opacities of lOssenkopf & Hennind )l994h by a factor of 1.5. 
This reduced value for the opacities has not been applied when deterniining 
the clump masses presented here and therefore we have rescaled the value 
given bvlKautfmann et al. 1 2010b ) of 870 by this factor to the value of 580 
given m Eqn. 5 (cf. lDunham et all201 id) . 
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Figure IS. The mass-size relationship of the ATLASGAL-MMB associated 
clumps. The yellow shaded region shows the part of the parameter space 
found to be devoid of massive star formation that satisfies th e relationship 
m(r) < 580 Mo (fieff/pc)' " (cf IKauffmann & Pillail [20101) .^ The green 
shaded region indicates the region of parameter space wh ere the young 
massi ve cluster progenitors are expected to be found (i.e., iBressert et al.l 
I201I . The dashed blue line shows the result of a linear least-squares fit to 
the resolved dust sources. The grey dashed line shows the sensitivity of the 
ATLASGAL survey and the upper and lower solid red line shows the sur- 
face densities of 1 gcm"^ and 0.05 gcm~^, respectively. The diagonal pink 
band fills the gas surface density (£(gas)) p arameter space betwee n 116 - 
129 Mo pc"- suggested bv lLada et al] <2010l) and Heiderman et afj j2010h . 
respectively, to be the threshold for "efficient" star fonnation. 



the less massive clumps (i.e., less than 1,000 M©) also have the po- 
tential to form massive stars, as long as they are relatively compact. 

Only 6 ATLASGAL-MMB associated sources are in the part 
of the parameter space that was found to be devoid of massive stars. 
All of these sources have been placed at the near distance and if the 
wrong distance has been assigned then this would explain their lo- 
cat ion in the mass-radius plot. Checki ng the confidence flag given 
bv iGreen & McClure-Griffith^ (1201 lb we find that three of these 
sources have been given a flag of b, indicating their distance assign- 



mass-radius requirement for mas- 



ments are less reliable, and may explain why these sources fail to 
satisfv lKauffmann et al.l (|20 
sive star formation. 

Turning our attention back to Fig.[T5] we see that the 
ATLASGAL-MMB data form a fairly continuous distribution over 
almost four orders of magnitude in mass and two orders of mag- 
nitude in radius. The dashed blue line overlaid on this plot shows 
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the result of a linear least-squares fit to the data. This fit provides 
a good description of the mass-size relationship for these objects 
(Log(Mdi,n,p) = 3.4 ± 0.013 + (1.67 ± 0.036) x Log(/?eff))- The up- 
per and lower red diagonal lines indicate constant surface densities, 
l(gas), of 1 gcm~^ and 0.05 gcm"^, respectively. These two lines 
provide fairly reliable empirical upper and lower bounds for the 
clump surface densities required for massive star formation. Fur- 
thermore, the lower bound of 0.05 g cm"^ provides a better con- 
straint than Kauffmann et al. for the high-mass end of the distribu- 
tion (i.e., R^g > -0.5 pc or M,.ii,„,p > -500 Mq). 

The t hick pink line shows the threshold derived bv lLada et al.l 

( l2010h and lHeiderman et al.l ( |2010|) (1 16 and 129 Mq pc^^, respec- 
tively; hereafter LH threshold) for "efficient" star formation. Above 
this threshold the observed star formation in nearby molecular 
clouds (d < 500 pc) is linearly proportional to cloud mass and 
negligible below. This line corresponds to an extinction thresh- 
old of Ak - 0.9 mag or visual extinction, Av, of ^ 8 mag. As 
discussed in the previous paragraph, the lower mass-radius enve- 
lope of the ATLASGAL-MMB associations is well modeled by 
S(gas) = 0.05 gcm"'^, which is approximately twice the LH thresh- 
old. However, we note that the star formation associated with the 
nearby molecular clouds used to determine the LH threshold is 
likely to be predominantly low-mass, while the threshold deter- 
mined from the ATLASGAL-MMB associations may be a require- 
ment for efficient formation of intermediate- and high-mass stars. 
This resul t suggests that a volume density threshold may apply 
(e.g.. |Parmentier et al.,2Ql ij) , however, with the data presented here 
it is not possible to determine the density distribution or the volume 
density at the scales involved with star formation. 

Given that this sample includes a large range of physical sizes 
(cores, clumps and clouds) and evolutionary stages (HMC, HMPO 
and UC H ii regions) it is somewhat suiprising to find such a strong 
correlation. However, as mentioned in Sections 4.2 and 4.3, the 
clumps appear to have a scale-free envelope structure that is not 
significantly changed as the embedded YSOs evolve towards the 
main sequence. 



5.1.2 Precursors to young massive clusters 

The green shaded area in the upper left part of Fig.|15 1 indicates 
the region of the mass-radius parameter space in which young 
mas sive proto-cluster (M PC) candidates are thought to be found 
(see iBressert et al.] 1201 2| for details). MFCs are massive clumps 
with sufficient mass that, assuming a fairly typical SFE (i.e., 
-30 per cent), have the potential to be the progenitors of future 
young massive clusters (YMCs) that have masses of > 10"^ Mq suc h 
as the Arches and Quintuplet clusters JPortegies Zwart et al.l201oh . 
We consider both starless and star-forming clumps to be MFCs as 
anything with a gas e nvelope is a possible p recursor. In an anal- 
ysis of the BGFS data. lGinsburg et al] l l2012h identified only three 
MFC candidates in a longitude range of ^ = 6-90° f rom a sample of 
-6,00 dust clumps. Using their detection statistics. lGinsburg et alj 
( l2012h go on to estimate the number of MFCs in the Galaxy to be 
< 10 ± 6, however, they state that a similar study of the southern 
Galactic plane is needed. 

Given that the lifeti me of the starless massive clumps is 
relati vely short (-0.5 Myr; iTackenberg et al.ll2012l : iGinsburg et al.l 
1201 2l) it is probably safe to assume that the ATLASGAL-MMB 
sample probably includes all of the likely MFC ca ndidates found 
in the 20° > ^ > 280° of the Galaxy. Applying the IBressert et alj 
( I2OI2I) criteria to our ATLASGAL-MMB sample, we have identi- 
fied 7 young massive proto-cluster candidates, one of which (i.e.. 
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Figure 16. Galactic latitude distribution of ATLASGAL sources (grey filled 
histogram) and ATLASGAL-MMB associated sources (blue histogram). 
The bin size used is 0.1°. 



AGAL 010.472-h00.027) is loc ated in the overlapping regions anal- 
ysed bv lGinsburg et al.l ( l2012h who also identified this source as an 
MFC candidate. The ATLASGAL clump names and their derived 
parameters a re tabulated in Table[3 ]alon g with the MFC candida tes 
identified bv lGinsburg et al.l (1201 2h and Longmore et al ] (l2012bl) . 

We note that the detection rate of young massive proto-cluster 
candidates identified from the ATLASGAL-MMB sample is twice 
that identified from the BGFS l lGinsburg et al. 2012). Both stud- 
ies covered a similar area of the Galactic disk (approximately 
30 per cent each between galactocentric radii of 1-15 kpc), how- 
ever, ATLASGAL covers a much broader range of Galactic lati- 
tude (1*1 < 1.5° and \b\ < 0.5° for ATLASGAL and the BGFS, re- 
spectively) and we note that three of the MFC candidates identified 
here have latitudes greater than 0.5°. The difference in MFC can- 
didate detection rates between the two surveys is likely a result of 
the difference i n latitu de coverage. It is therefore also likely that 
tai nsburg et al.l ( |20I2|) have underestimated the number of MFC 
candidates in the first quadrant by a factor of two. Combining the 
number of MFCs identified and the volume coverage we estimate 
the total number of MFC candidates in the Galaxy is < 20 ± 6. 

The number of MFC candidates is a factor of two larger 
than the number of YMCs known in the Galaxy, but similar to 
the number of embedded clusters in the Galaxy estimated by 
iLongmore et al.l f2012bl) (i.e., -25). However, this is almost an or- 
der of magnitude greater than the number of MFCs expected as- 
suming a YMC lifetime of 10 Myr and a formation time of - 1 Myr. 
This would suggest that either only a small number of the MFC can- 
didates will evolve into YMCs, which would require a SFE lower 
than the assumed 30 per cent, or there are many more YMCs that 
have as yet not been identified. 



5.2 Galactic distribution 

In the previous subsection we have shown that the clumps identi- 
fied from their association with methanol masers are likely to be 
involved with the formation of the next generation of massive stars. 
Massive star formation has been found to be almost exclusively 
associated with the spiral arms of nearby analogues of the Milky 
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Table 3. Derived parameters for massive proto-clusters. Masses have been estimated assuming a dust temperature of 20 K for all sources except 
G002. 5 3+00. 016 for which a range of temperatures between 19 and 27 K estimated from SED fits to individual pixels across the clump (see lLongmore et al] 
l2012bl for details). 



ATLASGAL Name 


Distance 


Effective Radius 


Clump Mass 


Reference 




(kpc) 


(pc) 


(Log(Mo)) 




AGALO 10.472+00.027 


11.0 


3.0 


4.77 




AGAL328.236-00.547 


11.4 


4.9 


4.99 




AGAL329.029-00.206 


11.7 


3.2 


4.68 




AGAL345.504+00.347 


10.8 


5.7 


4.93 




AGAL350.il 1+00.089 


11.4 


2.1 


4.55 




AGAL35 1.774-00.537 


17.4 


4.8 


5.43 




AGAL352.622-01.077 


19.4 


3.3 


4.79 




GOOO.253+00.016 


8.4 


2.8 


5.10 


2 


GO 10.472+00.026 


10.8 


2.1 


4.58 


3 


G043. 169+00.009 


11.4 


2.2 


5.08 


3 


G049.489-00.370'' 


5.4 


1.6 


4.68 


3 


G049.489-00.386'' 


5.4 


1.6 


4.72 


3 



References: (1) this work, (2) lLongmore et al.l l l2012"bl) . O^ lGinsburg et al.l l l2012h 

Notes: "These two sources are considered as a single MPC in the discussion presented bv lGinsburg et al.M2012l) . 
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Figure 17. Galactic longitude distribution of ATLASGAL sources (grey filled histogram) and ATLASGAL-MMB associated sources (blue histogram) are 
presented in the lower panel, while in the upper panel we present the fraction of ATLASGAL sources associated with a methanol maser The bin size used is 
2°. 



Way jKennicutl|[200^) and therefore the Galactic distribution of 
this sample of ATLASGAL-MMB associated clumps may provide 
some insight into where in the Galaxy massive star formation is 
taking place and its relation to the spiral arms. 

Figs.|16|and|17| we show the distribution in Galactic latitude 
and longitude of ATLASGAL sources and ATLASGAL-MMB as- 



sociations in our sample, in the range 20° > f > 280°. We find 
no significant difference between the ATLASGAL-MMB sample's 
latitude distributions and the overall ATLASGAL s ource distribu- 
tion; w hic h has been prev i ously commented on bv iBeuther et alj 
( I2OI2I) and lContreras et all ( l2013h and so will not be discussed fur- 
ther here. The distribution of source counts at 2-degree resolution 
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is shown in the lower panel of Figs.[T7] and the fraction of AT- 
LASGAL sources with MMB associations in the upper panel of 
this figure. While the total number of sources is dependent on both 
the local source density and the line-of-sight distribution, the ratio 
is independent of these and shows the incidence rate of methanol 
masers in dense clumps. 

The distribution of ATLASGAL sources as a function of 
Galactic longitude in Fig.[T7] (lower panel) reflects Galactic struc- 
ture, with a strong, narrow peak towards the Galactic Centre and 
a broader maximum centred around i ~ 330, which corresponds 
to a group of sources in the Scutum-Centaurus arm and possibly 
along th e Norma arm tangen t, which is close to the s ame direc- 
tion (see lBeuther et alj|2012l and IContreras et al]|2013l for further 
discussion). Since methanol ma sers are associated with the early 
stages of massive star formation dMinier et al.ll2003h . maser source 
counts can be taken as a measure of the massive star-formation rate 
(SFR). If the total maser counts trace the massive SFR on this spa- 
tial scale, then this fraction traces an analogue of the star-formation 
efficiency (SFE), i.e. the rate at which dense, submillimetre-traced 
clumps are producing massive YSOs, within the timescale appro- 
priate to the metha nol-maser stage of evolution (~2.5-4.5 x 10* yr; 
Ivan der Waldl2005h . The MMB fraction dips significantly at longi- 
tudes \t\ <4°, and perhaps within 10°, suggesting that the SFE 
is significantly reduced near the Galactic centre. Outside this nar- 
row zone, there is little evidence of any significant changes in SFE 
on these scales, associated with other features of Galactic struc- 
ture within the survey area. In particular, the ratio is more or less 
constant across both the peak in source counts seen at ^ ^ 330° 
mentioned above, and the clear drop in counts at larg er longitudes. 
Similar results were reported by iBeuther et al.l (|2012|) from a com- 
parison of GLIMPSE source co unts to the ATLASGAL longitude 
distribution. iMoore et al.l l l2012l) report increases in the large-scale 
SFE a ssociated with some spiral-arm structures, while Eden et alj 
( l2012h found no significant variations associated with major fea- 
tures of Galactic structure in the fraction of molecular cloud mass 
in the form of dense clumps. 

The lower SFE towards the Galactic centre inferred from the 
ratio of ATLASGAL-M MB associated clumps is supported by a 
recent study presented bv lLongmore et all ( l2012al) . These authors 
used the integrated em ission of the inver s ion transition of am mo- 
nia (1,1) from HOPS jWalsh et alJbol ll : IPurcell et al.ll2012h and 
the 70-500 /jm data from the Hi-GAL survey jMolinari et alj2010l) 
to trace the Galactic distribution of dense gas and compared this 
to star formation tracers such as the water, methanol masers and 
UGH n regions identified by HOPS, MMB a nd the Green Ban k 
Telescope HII region Discovery Survey (HRDs JBania et alj2010h . 
respectively, to derive the star formation rate per unit gas mass. 
The overall distribution of the integrated ammonia emission is very 
similar to that of the ATLASGAL source distributi on showing a 
very s trong peak towards the Galactic cent re (cf. Fig. 3 |Purcell et alj 
l2012h . Although lLongmore et al.l ( l2012ah found the ammonia and 
70-500 yum dust emission to be highly concentrated towards the 
Galactic centre they found the distribution of masers and H ii re- 
gions to_be_jBlatiyely uniform across the Galaxy. Measuring the 
SFR. lLongmore et al.l ( 1201 2al) found it to be an order of magnitude 
lower than what would be expected given the region's surface and 
volume densities. 

We are unable to independently estimate the absolute SFR 
from the MMB and ATLASGAL source counts. However, we can 
state that we find the fraction of clumps associated with the forma- 
tion of massive stars is a factor of 3-4 lower in the Galactic centre 
compared with that found for the rest of the Galaxy. The lower 



SFE seen towards the Galactic centre is likely a result of the ex- 
treme environment al conditions found in this region of the Galaxy 
dimmer etal .120 12h . 

In Fig.[T8] we present the 2-D Galactic distribution of 
ATLASGAL-MMB associations. In this figure we only include 
ATLASGAL-MMB associations that have masses and maser lumi- 
nosities for which the sample is complete across the Galaxy (i.e., 
M > 1,000 Mo and Lmmb ^ 1,000 Jy km s"' kpc^). We find that the 
positions of the ATLASGAL-MMB associations are in reasonable 
agreement with the main structural features of the Galaxy, as deter- 
mined by other tracers (taking into account the kinemat ic distance 
uncertainty due to peculiar motions is of order ±1 kpc: iReid et alj 
l2009l) . However, we also note that the degree of correlation is not 
uniform, with high source densities coincident with the near section 
of the Scutum-Centaurus arm and the southern end of the Galactic 
Bar. There is also a smaller number of sources that are distributed 
along the length of the Sagittarius arm. However, one of the most 
interesting and surprising features of Figs.[T7]and[T8]is the lack of 
any clear enhancement in ATLASGAL sources or ATLASGAL- 
MMB associations corresponding to the Scutum-Centaurus arm 
tangent at ^ ~ 3 10° -320°. The low number of clumps with or with- 
out masers associated with the southern Scutum-Centaurus am im- 
plies that this arm is not strongly forming stars of any type, and 
not just a lack of high-mass star formation. This is more surprising 
when we consider the distribution of CO emission, which is found 
to peak between { = 310 - 312° (see Figs. 6 and 7 presented in 
[Bronfman et aLlfTQSSl) . This would suggest that this arm is associ- 
ated with significant amounts of molecular material, but that the 
clump forming efficiency is relatively low for some reason. 

Another interesting feature seen in Fig. [18] is the relatively 
low number of ATLASGAL-MMB asscociated clumps located at 
large Galactic radii. Since we have only plotted clumps with masses 
above 1,000 Mq this low number of sources found outside the solar 
circle is not due to a lack of sensitivity. In Fig.[T9]we present the 
ATLASGAL-MMB source surface density as a function of Galac- 
tocentric radius. The three peaks in the Galactic radial distribution 
at ~3, 5 and 10 kpc correspond to the far 3-kpc spiral arm, a com- 
bination of the near section of the Scutum-Centaurus arm and the 
southern end of the Galactic Bar, and the Sagittarius spiral arm, 
respectively. There is a strong enhancement in the massive star- 
formation rate in the far 3-kpc arm with an otherwise relatively 
constant surface density between 3 and 7 kpc, however, the surface 
density drops significantly at large radii. 

There is clearly a significant difference between the surface 
density inside and outside of the solar circle. The co-rotation radius 
is approximately the same as the solar circle (i.e, Rqc = 8.5 kpc), 
which is where the spiral arms — in principle — are less important 
because the ISM is no longer being shocked as it runs into an arm. 
So one explanation could be that the spiral arms play an important 
role in creating the conditions required for massive star formation 
within the inner Galaxy (e.g., efficiently forming molecular clouds 
from the gas entering the spiral arm) compared to the outer Galaxy. 
The co-rotation radiu s is also where the m etallicity has been found 
to drop sharply (e.g.. iLepine et al.ll201 ll) and therefore metallicity 
may also play a role. 

Since the ATLASGAL-MMB sample is primarily tracing mas- 
sive star formation it is unclear whether this difference in the sur- 
face density of star formation between the inner and outer Galaxy 
is restricted to the massive SFR or if it is also found in the low- 
and intermediate-mass star formation; if it is a feature of the mas- 
sive SFR then it could have implications for the IMF. However, it is 
also possible that the galactocentric massive star formation surface 
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Figure 18. Galactic distribution of the ATLAS GAL-MMB sources with known distances and clump masses ^ 10^ Moand Lmmb > 1,000 Jy km s"' kpc^). We 
show the kinematic positions of our sample as red-in-black (yellow-in-black for the PMC candidates) circles, the sizes of which give an indication of their 
respective dust mass. In the upper right comer we give the masses for a sample of clumps. We have superimposed the ATLASGAL-MMB source distribution 
over a sketch of how the Galaxy is thought to appear if viewed face-on. This image has been produced by Robert Hurt of the Spitzer Science Center in 
consultation with Robert Benjamin and attempts to synthesise many of the key elements of Galactic structure using the best data currently available (Courtesy 
NASA/JPL-Caltech; see text for more details). The position of the Sun is shown by the small circle above the image centre. The Roman numerals in the comers 
refer to the Galactic quadrants and the two white solid lines originating from the location of the Sun enclose the region of the G alactic Plane overlapped by 
the MM B and ATLASGAL surveys, while the two white dashed lines indicate the region towards the Galactic centre excluded bv lGreen & McClure-Griffithsl 
i201lh as kinematic distance in this region are unreliable. The dot-dashed circles represent the locus of tangent points and the solar circle. 



density is simply reflecting the underlying distribution of molecu- 
lar and atomic gas in the Galaxy. This is something we will revisit 
when distances and masses are available for a larger fraction of the 
ATLASGAL compact source catalogue. 

The surface density distribution presented in Fig.[T9]has been 
determined using only the ATLASGAL-MMB associations that are 
above the clumps mass and maser luminosity completeness limit. 
By multiplying the surface density in each of the bins in this plot 
by the area of the bin annulas we estimate the total number of mas- 
sive star forming clumps with masses > 1000 Mq in the Galaxy to 



be ~560. The sample presented here therefore represents approxi- 
mately 50 per cent of the whole Galactic population of massive star 
forming clumps that are associated with a methanol maser. 

We can take this analysis a step further and estimate the con- 
tribution these star forming clumps make to the Galactic SFR. If we 
take the median clump mass of ~3000 Mq and a SFE of 30 per cent 
then combined we estimate these clumps will produce clusters with 
a total stellar mass of ~ 5 x 10^ Mq. The formati on time for the 
most m assive stars in these clus ter was calcul a ted bylMolinari et"^ 
( I2OO8I) to be -1.5x10^ yr and lOavies et alj ( I2OIII) estimated the 
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Figure 19. Surface density of ATLASGAL-MMB associated sources as 
a function of galactocentric distance. As witli Fig.[T8] we only include 
ATLASGAL-MMB associations that have masses and maser luminosities 
for which the sample is complete across the Galaxy (i.e., M > 1,000 Mq 
and Lmmb ^ 1 ,000 Jy km s" ' kpc^ ). The bin size is 0.5 kpc. 

combined statistical lifetime of the MYSO and UCH ii region 
stages to be several 10' yr, and therefore a reasonable lower limit 
for the cluster formation time is likely to be ~0.5 Myr. However, 
most of the total cluster stellar mass comes from lower-mass stars, 
which take longer to form and will increase the cluster formation 
time to perhaps 1 Myr. Using these two cluster formation times as 
an upper and lower limit we estimate that these methanol maser 
associated clumps have combined SFR of between 0.5-1 M© yr"', 
and therefore could be responsible for u p to half the current Galac- 
tic SFR ~2Mo yr ' l lDavies etalhouh . 

5.3 Luminosity-mass correlations 

In Sect. |4.5| we found a weak correlation between the clump mass 
and maser luminosity and speculated that this may be related to 
higher luminosities of the embedded young stars. In this section 

we will explore this pos sible link. 

In a recent paper iGallawav et al.l | |2013|) reported the posi- 
tional correlation of M MBs with other samp les of young high mass 
stars such as the RMS jUrquhart et al.l2008h and with EGOs which 
are thought to trace shocked g as associated with the outflows of 
MYSOs (e.g., Cv ganowski et a l. 2008, 20(3)- They identified 82 
RMS sources (UC H it regions and/or YSOs) within 2" of an MMB 
source and found that all of these have luminosities consistent with 
the presence of an embedded high-mass star. We have repeated the 
spatial correlation described by Gallaway et al., however, we re- 
laxed the matching radius to 20" in order to determine the standard 
deviation of the offsets and set the association radius accordingly; 
this identifies 137 possible matches. In Fig.|20]we present a plot of 
the surface density as a function of angular separation between the 
matched RMS and MMB sources. The distribution is peaked at off- 
sets between and 1" and rapidly falls off to approximately zero 
for separations greater than 4". We consider ATLASGAL-MMB 
sources to be associated with an RMS source if the offset is less 
than 4"; this identifies 102 reliable matches, which corresponds to 
approximately 20 percent of the MMB sample in the 20° > > 10° 
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Figure 20. Surface density of RMS-MMB associations as a function of 
separation. We have truncated the j:-axis of this plot at 15". The bin size is 
1". 
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Figure 21. Bolometric luminosity distribution of the RMS sources associ- 
ated with a methanol maser. The bin size is 0.5 dex. 

and 350° > C > 280° region (the RMS survey excluded sources 
within 10° of the Galactic centre due to the increased background 
confusion and larger distance uncertainties). 

Bo lometric luminosities have been estimated for the RMS 
sources jMottram et al.l2010l l2oT I]) and thus allow us to investigate 
the luminosity distribution of these RMS-MMB associations. The 
fluxes used to estimate these luminosities are effectively clump- 
average values and are therefore a measure of the total cluster lumi- 
nosity. However, if we assume that the luminosity is emitted from 
zero age main sequence (ZAMS) stars then the cluster luminosity 
is dominated by the most massive star in the cluster and so can be 
used to probe the correlation between it and the methanol maser 
Since many of these masers are associated with UC H it regions the 
ZAMS assumption is reasonable. 

The luminosities range from ~ 100 to 10^ Lq with a peak in the 
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Figure 22. The bolometric luminosity plotted against clump mass (upper 
panel) and methanol maser luminosity (lower panel). Sources classified as 
YSOs and UCHII regions are shown as filled and open circles, respectively. 
The red line show the result of a lineal' least-squares log-log fit to the whole 
sample. 



distribution at 10"""^" Lg (see Fig.l21b. The RMS survey is com- 
plete for emb edded young mass ive stars with luminosities a few 
times lO'' Lg jPavies etalll201lh and so the fall off in the distri- 
bution for luminosities lower than this value is due to incomplete- 
ness. The minimum luminosity of a high-r nass star is ~ 1,000 Lm 
(corresponding to a star o f spectral type B3: ISoelmi-Vitens^ l 19811 : 
iMevnet & Maedei|[2003l) and so the fact that most of the RMS- 
MMB sample (~90 per cent) have luminosities above this supports 
their association with high-mass star formation. 

While the majority of the MMB sources appear to be associ- 
ated with massive stars we do find a handful of sources with lumi- 
nosities towards the lower end of the distribution (i.e., 100-500 Lq). 
EiTors in the kinematic distance solution could possibly account 
for a few of these but it is unlikely to explain them all. It there- 
fore seems likely that these methanol masers are associated with 
intermediate-mass stars with masses of several sol ar masses. This is 
consistent with the findings o f lMinieretal1 ( l2003h who were able to 



put a lower limit of ~3 Mq for strong maser emission. It is possible 
that these are intermediate mass protostars that are still in the pro- 
cess of accreting mass on their way to becoming a high-mass proto- 
star. This is consistent with the mass-size relationship discussed in 
the previous subsection, however, as we will show in the next few 
paragraphs, there is a strong correlation between clump mass and 
bolometric luminosity and therefore these low luminosity sources 
also tend to be associated with the lower mass clumps and would 
require a high SFE (>50 per cent) to form a massive star. 

In the upper and lower panels of Fig.|22]we plot the bolometric 
luminosity as a function of clump mass and methanol maser lumi- 
nosity for these ATLASGAL-MMB-RMS associated sources, re- 
spectively. In these plots we indicate the evolutionary type as clas- 
sified by the RMS team, however, a KS test is unable to reject the 
null hypothesis that the two source types are drawn from the same 
populations and therefore we do not distinguish between them here. 
Again calculating the partial coefficients to remove the dependence 
of these parameters on distance we find a strong correlation be- 
tween the bolometric luminosity and clump mass (r=0.78 and p- 
value « 0.01) and a weaker correlation between the bolometric and 
maser luminosities (r=0.42 and p- value 0.01). The linear log- 
log fit to the bolometric luminosity and clump mass gives a slope 
of 0.94 ± 0.04 and so not only are these two parameters strongly 
correlated but their relationship is very close to linea r. This value is 
slightl y shallower than the slope of ~ 1.3 reported bv lMolinari et al.l 
|2008 !) from their analysis of a sample of 42 regions of massive star 
formation. However, we note that they use IRAS fluxes to deter- 
mine source luminosities and this will tend to overestimate fluxes 
for the more distant sources, which are typically also the most mas- 
sive sources. This would lead to a steepening of the luminosity- 
mass slope when compared to the RMS luminosities, which used 
the higher resolution 70 /jm MIPSGAL band to compute luminosi- 
ties and therefore should be more reliable. 

Assuming that the clump mass remains relatively constant 
through the early embedded stages of massive star formation while 
the luminosity increases as both the accretion rate and the mass of 
the protostar increase, the source should move vertically upwards in 
the mass-luminosity plot. Once the H ii region has formed and starts 
to disrupt its host clump the source's luminosity will be constant 
(luminosity of the ZAMS star ionizing the H ii region) and its mass 
will begin to decrease resulting in it moving horizontally to the 
left. Since there is no significant difference between the luminosity- 
mass distributions of MYSO and UC H ii region stages we can as- 
sume they have similar ages; at both of these stages core hydrogen 
burning has begun, however, the UC H n regions are slightly more 
evolved and have started to ionise their surroundings. The fit to the 
data (solid red line in the upper panel of Fig.l22t is a crude ap- 
proximation of the transition between protostellar and H ii region 
evolution. 

As Gallaway et al. point out, the MMB sources are likely to 
be associated with a broad range of evolutionary states from HMCs 
to UC H II regions with their luminosity increasing through accre- 
tion as they develop. As previously mentioned only ~20 per cent 
of ATLASGAL-MMB associated sources are associated with the 
later stages (i.e., YSO and UCH ii) and therefore the majority are 
younger and less evolved embedded objects whose luminosity is 
not yet sufficient to dominate that of its associated cluster. All of 
these sources will be located below the red line with more evolved 
H II regions dominating the upper left part of the parameter space. 

We have found the mass and bolometric luminosity to be 
weakly correlated with the methanol maser luminosity. The mass 
and bolometric luminosity are strongly correlated with each other 
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Figure 23. Examples of the local mid-infrared environment found towards the unassociated methanol masers. Images and symbols are as described in the 
caption of Fig.[T] The contour levels stait at lo" and increase in steps of Icr to emphasis any weak compact or diffuse emission in the vicinity of the methanol 
maser. In the upper panels we present images of MMB sources that show extended 8/jm emission commonly associated with star forming regions, while in 
the lower panels we present images of MMB sources that appear to be isolated and therefore possible evolved stars. 



(i.e., r = 0.78) and have very similar correlation values when com- 
pared to the maser luminosity (i.e., r = 0.42 and 0.44, respectively), 
which might make the underlying cause of the correlation hard to 
distinguish. However, since the maser arises from a population in- 
version pumped by far-infrared and submilhmetre photons emitted 
from the embedded YSO, it is the intrinsic luminosity of the YSO 
that is more likely the fundamental cause of these correlations. The 
weaker correlation between the bolometric and maser luminosities 
is probably related to the fact that the former is a measure of the 
total luminosity of the embedded proto-cluster, whereas the latter 
is more likely to be driven by a single cluster member, and no doubt 
maser variability also plays a role. 

A similar correlation between b olometric luminosity an d wa- 
ter maser luminosity was reported by lUrquhart et al.l fcOl Ibl) from 
a single dish survey of MYSOs and UC H ii regions identified by 
the RMS survey. Having found a correlation between the methanol 
and water maser luminosities for these two advanced stages we 
would conclude that the strength of the maser emission is dom- 
inated by the energy output of the central source and not driven 
by source evolution as suggested by Breen et al. (see discussion 
in Sect. |4.4.rt . One caveat to this is that both models of massive 
star formation (i.e., monolithic collapse and competitive accretion) 
predict that massive stars gain their mass through accretion from a 
circumstellar disk, and therefore their maser luminosity is likely to 
increase proportionally with their mass and luminosity. However, 
the large range of possible final masses of OB stars renders the 
maser luminosity insensitive to the current evolutionary stage of a 
particular stars in the same way that the bolometric luminosity is 
insensitive to evolutionary state of the embedded object. 



5.4 The nature of the unassociated MMB sources 

In Sect. l3.3l we identified 43 methanol masers that were not 
matched to an ATLASGAL source. In this section we will exam- 
ine the available evidence to try and investigate the nature of these 
sources. 

We begin by inspecting the ATLASGAL maps at the locations 
of the 43 MMB sources. In many cases, the position of the maser is 
coincident with either low-surface brightness, diffuse submillime- 
tre emission, which would have been filtered out by the background 
subtraction used in the source-extraction process, or weak compact 
emission that fell below the detection threshold used in the source 
extraction. We have measured the peak 870 yum flux at the position 
of the MMB source and estimated the 3(T noise from the standard 
deviation of nearby emission-free regions in the map, and have used 
the higher of these two values as the upper limit to the submillime- 
tre flux. We have used these values with their assigned distances (as 
discussed in Sect. l4.U to estimate an upper limit for their masses 
using Eqn. 1 and again assuming a dust temperature of 20 K. These 
results are summarised in Table|4] 

As previously mentioned, it is widely accepted that methanol 
masers are almost exclusively associated with high-mass star- 
forming regions. This is supported by many of the findings pre- 
sented in this paper. One explanations is that these unassociated 
MMB sources are located at larger distances and that their sub- 
millimetre emission simply falls below the ATLASGAL detection 
sensitivity. Indeed, the median distance for the unassociated maser 
sample is ~12kpc, which is much larger than the median value of 
~5 kpc found for the ATLASGAL-MMB associations and the KS 
test showed the two distributions to be significantly different (see 
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Table 4. MMB sources without a counterpart in the ATLASGAL catalogue. 
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Fig.|7]for comparison of distance distributions and Sect. l4. ll for dis- 
cussion). Looking at the estimated upper limits for the mass we 
find that they are all significantly lower than the 1 ,000 M© assumed 
to be the minimum required for massive star formation. However, 
as shown in Sect. |5.1| it is at least feasible for less massive clumps 
to form massive stars, but in the majority of cases these sources 
would need to be significantly smaller than the beam. The only 
caveat is that for all of the other ATLASGAL-MMB associations 
we found a scale-free envelope, but these sources would need to 
be relatively discrete and isolated clumps that are not embedded 
in a larger structure. Alternatively, it is also possible that these 
methanol masers are associated with embedded sources that will 
go on to form intermediate-mass stars, in which case the low val- 



ues obtained for the upper limits to the masses may not be all that 
important. 

Assuming this hypothesis is correct then, even though the dust 
emission is too weak to be detected, we might expect to see asso- 
ciated mid-infrared emission from diffuse nebulosity and evidence 
of extinction from dark lanes of dust often seen towards sites of 
star formation, bf Approximately 80 per cent of these sources are 
associated with mid-infrared emission (these are indicated by su- 
perscripts given by the MMB name in Table[4) , which is similar to 
the proportion found bv lGallawav et al.ll2013l for the whole MMB 
catalogue (i.e., 83 per cent). We present a sample of these false- 
colour mid-infrared images in Fig.|23]again created by combining 
data extracted from the GLIMPSE archive. 



The upper middle panel of this figure shows the mid-infrared 
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image of the MMB source MMB303.507-00.721, which is located 
at a distance of ~llkpc. This source is found to have extended 
8 /jm emission commonly associated with star-forming regions and 
there is evidence of compact dust emission to the south-west of 
the mid-infrared emission. From a visual inspection of the mid- 
infrared emission we estimate that in approximately a third of cases 
the emission is consistent with this maser being associated with a 
more distant star-formation site, however, this is probably a lower 
limit. 

So for a significant fraction of these unassociated methanol 
masers it is likely that the lack of dust emission can be explained by 
them being more distant and their associated dust emission falling 
below the ATLASGAL detection threshold. However, there is also 
a small number of sources for which this explanation is not satis- 
factory. There are three sources that are located at relatively near 
distances (i.e., MMB292.468-h00.168, MMB3 11. 55 1-00.055 and 
MMB3 12.702-00.087) where we would expect to have detected 
their dust emission. There are two more that are located at the far 
distance, but where the far distance allocation puts them much far- 
ther from the Galactic mid-plane than expected for star forming re- 
gions (i.e., MMB3 11.729-00.735 and MMB33 1.900-01. 186 that 
have z distance of - 1 83.6 and -244.5 pc, respectively), which casts 
some doubt on the distance allocation if indeed these are star form- 
ing. However, there is another intriguing possibility which is that 
these masers may arise in the circum stellar shells associated with 
evolved stars (e.g.. lWalsh et al.f2003l) . 

A search of the SIMBAD database revealed that only ten of 
these methanol masers were previously known: one identified in 
the IRAS point-source catalogue; two detected in the BGPS, and so 
it is likely that their dust emission falls below our detection thresh- 
old but is detected by the BGPS due to the ir superior low surface- 
brightness sensitivity; and 7 are included in lRobitaille et al.l l l2008h 
intrinsically red GLIMPSE source catalogue, six of which they 
classify as YSOs, and one (MMB328.385-l-00.131) they classified 
as a possible asymptotic giant branch star. In the lower right panel 
of Fig.|23] we present the mid-infrared image of the MMB source 
MMB328. 385+00. 131 that shows only a single point source coin- 
cident with the position of the methanol maser. This point source is 
isolated in the image and the lack of any extended 8 yum emission 
or extinction feature, that are commonly associated with star form- 
ing re gions, along with the classification made bv iRobitaille et al] 
( l2008h from its mid-infrared colours, would suggest this maser 
might be associated with an evolved star. However, this could sim- 
ply be due to a chance alignment along the same line of sight and so 
this association requires further investigation to test its reliability. 

Currently, there is not enough complementary data available 
for all of these sources to be able to properly evaluate these two 
explanation s. The Hi-GAL surv ey of the inner Galactic plane at 
70-500 yum l lMolinari et alj|2010l) will provide a way to definitively 
test these two possibilities (e.g.. lAnderson et^l20I2l) and will be 
discussed in a future publication. 



6 SUMMARY AND CONCLUSIONS 

The ATLASGAL survey JSchuUer et af] 120091; 280° < f < 60°) 
has identified the Galactic distribution of dust through its thermal 
emission at 870yum and is complete to all massive clumps above 
1,000 Mo to the far side of the inner Galaxy (~20kpc). In total the 
ATLASGAL survey has identified some 12,000 compact sources, 
many of which have the potential to form the next generation of 
massive stars. Methanol masers have been found to be associated 



with high mass star formation and we have therefore taken advan- 
tage of the availability of an unbiased catalo gue of these objects 
compiled by the methanol multibeam (MMB; ]Greedl2009l) survey 
team to identify a large sample of high-mass star forming clumps. 

Cross-matching these two surveys we have identified 577 
ATLASGAL-MMB associated clumps within the overlapping re- 
gion of both surveys (i.e, 280° < € < 20° and \b\ < 1.5°) with 
two or more methanol masers being detected towards 44 clumps. 
We find ~90 per cent of the matches are within 12" (~3cr) of the 
peak of the submillimetre emission revealing a strong correlation 
with column density and the location of a methanol maser within 
the clumps. We fail to identify any significant 870 //m emission to- 
wards 43 MMB sources. 

Ass uming a dust temperature of 20 K a nd using distances pro- 
vided bv lCreen & McClure-Griffithsl ( l20Ilh . and derived here, we 
are able to estimate the clump masses and radii, column and volume 
densities, and methanol maser luminosities for almost 500 of the 
ATLASGAL-MMB associations. We find we are complete across 
the Galaxy to all dust clumps with masses larger than 1,000 Mq 
hosting a methanol maser with luminosities > 1 ,000 Jy kpc^ . We use 
these parameters to investigate the link between methanol masers 
and massive star forming clumps and as a probe of Galactic struc- 
ture. Our main findings are as follows: 

(i) The clump radii cover a range from 0.1 to several parsecs 
with the larger clumps generally found at larger distances. With 
a median aspect ratio of 1.4 the ATLASGAL-MMB associations 
are fairly spherical centrally condensed structures, however, with a 
median K-factor of ~5 a significant amount of their mass is located 
outside the central region. The position of the methanol masers is 
strongly correlated with the peak column density at the centre of 
the clumps. We find no correlation between the aspect ratio and Y- 
factor with distance, which suggests that the envelope structures of 
these massive star forming clumps are scale-free. The apparently 
simple clump structure (with masers at the central col density peak 
and scale-free radial structure) suggests the formation of one cen- 
tral stellar cluster per clump. Formation of multiple clusters might 
be expected to be accompanied by more complex clump structure. 

(ii) We are complete to all dust clumps harbouring a methanol 
maser with masses over 1,000 M© across the inner Galactic disk 
(i.e., ~20kpc). The median clump dust mass (~3,000Mo) is sig- 
nificantly larger than the completeness level, which confirms that 
methanol masers are preferentially associated with massive clumps. 
Furthermore, assuming a Kroupa IMF and a star formation ef- 
ficiency of 30 per cent we find that 72 per cent of these clumps 
(i.e., Mciump > 1,000 Mq) are in the process of forming clusters 
hosting one or more 20 M© star(s). Although 28 per cent of the 
clumps have masses lower than 1,000 Mq we find these are also 
more compact objects (^0.3 pc) that are likely to form either sin- 
gle stars or small multiple systems of bound stars that are also very 
likely to include a massive star. This is supported by the empirical 
mass-r adius criterion for massive star formation dKauffmann et alj 
l20I0bl) that shows that 97 per cent of ATLASGAL-MMB associa- 
tions have masses and sizes that are consistent with all other known 
massive star forming clumps, including the lower mass clumps. We 
conclude that the vast majority of clumps associated with methanol 
masers are in the process of forming high-mass stars. 

(iii) Inspecting the mass-radius relation for the ATLASGAL- 
MMB associations we find that a surface density of 0.05 g cm"^ 
provides a better estimate of the lower envelope of the distribution 
for masses and radii great er than 500Maand 0.5 pc , respectively, 
than the criterion given bv lKauffmann et all JioiObl) . This surface 
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density threshold corresponds to a Aj^- ~ 2 mag or visual extinction, 
Av, of ^ 16 mag, which is approxima tely t wice the required t hresh- 
old determined bv lLada etid] ( l2010l) and iHeiderman et"al ](|20i3) 
for "efficient" low-mass star formation. This would suggest that 
there is a clear surface density threshold required for clumps be- 
fore star formation can begin but a higher threshold is required to 
for more massive star formation. 

(iv) Testing the evolu tionary trend reported in the literature (i.e. 
iBreenetalJboill |2012|) between the methanol maser luminosity 
and clump-averaged volume density we fail to find any correlation. 
Although we are able to reproduce the results of these previous 
studies we find that both parameters have a strong dependence on 
distance and that once this is removed the correlation between them 
drops to zero. However, we do find a the bolometric and methanol 
maser luminosities are correlated with each other. 

(v) We have identified seven clumps that have masses large 
enough to be classified as massive protocluster (MPC) candidates 
which are expected to forni the next generation of young massive 
clusters (YMCs) such as the present day Archers and Quintuplet 
clusters. Using the Galactic plane coverage of this study and the 
number of MPC candidates detected we estimate the Galactic pop- 
ulation to be <20±6. This value is twice as many as previously esti- 
mated and similar to the number of currently known YMCs, which 
would suggest that only a few of these MPC candidates will suc- 
cessfully convert the 30 per cent of their mass into stars required to 
form a YMC. 

(vi) The Galactic distribution reveals the Galactic centre region 
having a significantly lower star formation efficiency (SFE), ffian 
the rest of the Galaxy covered by this survey, which is broadly flat 
even towards the spiral arm tangents. The lower SFE is probably a 
reflection of the much more extreme environment found in the cen- 
tral region of the Milky Way. Interestingly we find no enhancement 
in either the ATLASGAL or ATLASGAL-MMB source counts in 
the direction of the Scutum-Centaums arm tangent, from which we 
conclude that this arm is not actively forming stars of any type. 

(vii) The galactocentric distribution reveals very significant dif- 
ferences between the surface density of the massive star formation 
rate between the inner and outer Galaxy. We briefly speculate on 
possible explanations, however, it is clear further work is required 
before the reasons behind this difference can be properly under- 
stood. Using the surface density distribution with the completeness 
levels applied we estimate the total Galactic population to be ~560, 
which means that the sample presented here represents approxi- 
mately 50 per cent of the whole population. We estimate the star 
formation associated with these methanol maser associated clumps 
may contribute up to 50 per cent of the Galactic star formation rate. 

(viii) Bolometric luminosities are available from the literature 
for ~100 clumps and these range between ~100 to 10* Lq with 
the distribution peaking at ~10^ Lq. This confirms the association 
between methanol masers and massive young stars for 90 per cent 
of this sample of clumps, but also reveals that there are some 
masers associated with intermediate-mass stars. For lower luminos- 
ity clumps (i.e., M^unip <1,000Lq) these may be intermediate-mass 
protostars that are still accreting mass and will eventually go on to 
form a high-mass star, however, these tend to be associated with the 
lower mass clumps and would therefore require a SFE >50 per cent. 
It therefore may be the case that a small number of these clumps 
are destined to only form intermediate-mass stars. 

(ix) We investigated the available evidence for the 43 methanol 
masers towards which no 870 /jm emission has been detected and 
concluded that while perhaps 50 per cent may be more distant star 
forming regions where the dust emission has simply fallen below 



the ATLASGAL surveys sensitivity, the nature of the other masers 
is yet to be determined. 

This is the first of a series of three papers planned to use the 
ATLASGAL survey to conduct a detailed and comprehensive in- 
vestigation of high-mass star formation. The main aim of these pa- 
pers is to use the unbiased nature of the dust emission mapped by 
ATLASGAL over the inner Galactic plane to connect the results 
derived from different high-mass star formation tracers. In subse- 
quent papers we will investigate the dust properties of an unbiased 
sample of ultra-compact HII regions identified from the CORNISH 
survey and a complete sample of massive YSOs identified by the 
RMS survey. 
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